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I. INTRODUCTION 



This is the fourth of a series of papers on the one dimensional noisy Burgers equation for 
the slope field of a growing interface. In paper I jl] we discussed as a prelude the noiseless 
Burgers equation P, [J in terms of its nonlinear soliton or shock wave excitations and per- 
formed a linear stability analysis of the superimposed diffusive mode spectrum. This analysis 
provided a heuristic picture of the damped transient pattern formation. As a continuation 
of previous work on the continuum limit of a spin representation of a solid-on-solid model 
for a growing interface jj], we a_p_plied in paper II [f| the Martin-Siggia-Rose formalism (| 
in its path integral formulation 3,0,0] to the noisy Burgers equation and discussed 

in the weak noise limit the growth morphology and scaling properties in terms of nonlinear 
soliton excitations with superimposed linear diffusive modes. In paper III ^| we pursued a 
canonical phase space approach based on the weak noise saddle point approximation to the 
Martin-Siggia-Rose functional or, alternatively the Freidlin-Wentzel symplectic approach 
to the Fokker-Planck equation |l3L Il4| . This method provides a dynamical system theory 
point of view to weak-noise stochastic processes and yields direct access to the probability 
distributions for the noisy Burgers equation; brief accounts of paper II and III appeared in 
Q and Q. 

Far from equilibrium phenomena are common, including turbulence, interface and growth 
problems, chemical reactions, and a host of other phenomena bordering on biology, sociology 
and economics. Unlike equilibrium phenomena the nonequilibrium cases are not very well 
understood and constitute a major challenge in modern statistical physics. Here the Burgers 
equation provides in many respects the simplest continuum description of a nonlinear initial 
value problem in the noiseless case and an open driven nonlinear system in the noisy case, 
exhibiting scaling and pattern formation. 

The noisy Burgers equation for the local slope u(x, t) of a growing interface analyzed in 
papers II and II has the form 

(J^ - XuV^j u = vV 2 u + V77; (1.1) 

here expressed as manifestly invariant under the slope-dependent nonlinear Galilei transfor- 
mation 

x — > x — Xu t, (1.2) 
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u — > u -\- uq, (1.3) 



and is equivalent to the much studied Kardar-Parisi-Zhang (KPZ) equation |17I . Il8j | for the 
height h(x,t), u = V/i (in a co- moving frame), 

^ = ^h + ^Vhf + r ] . (1.4) 

The growth equations ()1.1|) and (jl.4j) are driven by short-range correlated Gaussian white 
noise rj determined by the correlation function 

(r]r])(xt) = A8(x)6(t), (1.5) 

characterized by the noise strength A. In Eqs. (jl.l|) and (jl.4j) the damping constant 
or viscosity v measures the strength of the linear damping term, whereas A control the 
nonlinear growth or mode coupling term. 

From the analysis in papers II and III it follows that the stochastic nonequilibrium prob- 
lem determined by Eqs. (jl.lj) and (jl.5j) in the singular weak noise limit A — ► can be 
replaced by two coupled deterministic Galilean invariant mean field equations coupling the 
slope field u to the canonically conjugate noise field p 

(^-Awvjw = zA7 2 w-V 2 p, (1.6) 

\j^-\uv\p =-vV 2 p. (1.7) 

In the path integral analysis in paper II Eqs. (jl.fil - II. 7j) are saddle point equations for 
the extremal path in the classical limit A — > 0; in the canonical phase space approach in 
paper III they are classical canonical field equations determining orbits in an associated 
phase space. This doubling of dynamical variables in the deterministic description was also 
encountered in the spin model discussed in Q|. The noise variable rj in Eq. ()l.lj) emerges as 
the canonically conjugate momentum variable p coupling to u. 

As discussed in papers II and III, see also Q, the field equations (jl.fij) and (jl.7j) in 
addition to linear diffusive modes also support two distinct soliton modes or domain walls, 
in the static case of the kink-like form 

Ug(x) = fj J uta.nhk s (x — xq) , k s = \uj2v. (1.8) 

Here k s sets the inverse length scale, = ± is a parity index and the soliton is centered 
at Xq. For fi = +1 we have the right hand soliton which is also a solution of the noiseless 



Burgers equation for 77 = in Eq. or for p = in Eq. (|1.6|) . For /i = — 1 we obtain 

the noise-induced left hand soliton, a new solution of the coupled equations. The associated 
noise field is p s = for the noiseless // = +1 soliton; for the noisy soliton for // = —1 we 
have 



2i/«tanh 



;i.9) 



modulus a constant. 

In the noiseless Burgers equation the transient pattern formation is described by Galilee- 
boosted right hand solitons connected by ramp solutions with superimposed damped linear 



modes 
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2l| . In the noiseless KPZ equation for the height h this pattern corresponds 



to smoothed downward cusps connected by parabolic segments with superimposed linear 
modes However, in the noisy case the doubling of soliton solutions alters the 

morphology completely. Here the amplitude-matched Galilee-boosted right and left hand 
solitons provide a many body description of a stationary growing interface. On the soliton 
gas is superimposed a gas of linear modes which in the linear Edwards- Wilkinson case 2^| 
for A = becomes the diffusive modes of the noise-driven diffusion equation. 

The canonical phase space approach expounded in paper III moreover provides a de- 
terministic dynamical system theory description of a growing interface. With an orbit in 
canonical phase space from an initial configuration u l to a final configuration w f traversed in 
time T, determined as a specific initial-final value solution of the field equations and 
(11.71) . we thus associate the action 5* 



S{u\u\T) 



ul ' T , , ( du „/ 

atax p— ri 

uKo \ at . 
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More explicitly, the transition probability from an initial configuration u l to a final configu- 
ration u { in time T is determined by 



P(u\u\T) = n(T)- L exp 



-^5( M f ,<T) 



Q(T) = J n*/exp 



1.11) 

;i.i2) 



where we have introduced the dynamical partition function fi(T), arising from the normal- 
ization condition / T7 rfM f -P(w f , u\ T) = 1. Likewise, the stationary distribution is associated 
with an infinite-time orbit from u 1 to u l and is given by 



P st (u { )= lim P(u { ,u\T), 

T— >oo 



;i.i3) 



and for example the slope correlation (the second moment of P) in the stationary regime by 
the expression 



(uu)(x,T) = / Udu'du 1 u t (x)u 1 (0)P(u\u\T)P st (u 1 ). 



(1.14) 

The action is a central concept in the weak noise canonical phase space approach and provides 
a dynamical weight function and selection criteria for a dynamical nonequilibrium process 
in a similar manner as the energy E in the Boltzmann-Gibbs factor P oc exp[— [3E], (f3 is the 
inverse temperature) for equilibrium processes. The action moreover implies an underlying 
principle of least action and the Hamilton entering S, yielding the field equations ()1.6|) and 
(jl.7|) . is given by 



H = J dx p{uV 2 u + XuVu - (1/2) V 2 p), 



;i.i5) 



where H = f dx H. 

In addition to the conserved Hamiltonian or energy the translational invariance of H (as- 
suming periodic boundary conditions for u and for p, modulus a constant) implies conserva- 
tion of momentum II. Moreover, the conserved noise in Eq. (jl.lj) . corresponding to the term 
V 2 p in Eq. (jl.6j) . yields the local conservation law du/dt+Vj = 0, j = — v'Vu+'Vp— (X/2)u 2 , 
implying the conservation of the integrated slope field or height offset. The two additional 
conserved quantities are thus given by 



II = J dx uVp, 
M = dx u. 



(1.16) 
(1.17) 

We note that the integrated noise field P is not conserved for A ^ 0. According to Eq. (jl.6j) 
P evolves according to, see also 0, 



dP 
~dt 
P = 



The long time-large distance scaling properties of a growing interface is a 



An, 

dx p. 



issue which has been addressed extensively |Tj 
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(1.18) 
(1.19) 
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For the width w(L,t) of an interface of size L the dynamical scaling hypothesis 
asserts that w = L~^G(t/ L z ) which for the stationary slope correlations corresponds to the 



asymptotic scaling form 

(uu)(x, t) = x 2C ~ 2 F ) , (1.20) 

with roughness exponent ( = 1/2, dynamical exponent z, and universal scaling function 
F(w). In one dimension the scaling exponents for the noisy Burgers equation are known 
exactly Q, The roughness exponent ( = 1/2 follows from the known stationary distri- 
bution, an effective fluctuation-dissipation theorem 



— — / dx u 2 



P st (u) oc exp 

whereas the dynamic exponent z = 3/2 is a consequence of the scaling law 



1.21) 
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C + z = 2, (1.22) 



. It was an important result of the analysis in papers 



implied by Galilean invariance 

II and III, see also j^, that the dynamical exponent z = 3/2 also enters in the dispersion law 
E oc IF for the noise-induced left hand soliton and thus is a feature of the gapless nonlinear 
excitations providing the many body description of a growing interface. 

The description of the stochastic nonequilibrium dynamics of a growing interface can be 
accessed on two levels: The Langevin level defined by Eqs. (jl.lj) and (jl.5j) or the Fokker- 
Planck level (or Master equation level for discrete models) characterized by the Fokker- 
Planck equation associated with the Burgers equation, 

dP 

A— = -HP. (1.23) 

Here the Hamiltonian or Liouvillian H is given by Eq. (jl.lfrj) . with the momentum variable 
p interpreted as the functional derivative p = A5/Su, see also |4(. 

On the Langevin level the growth problem is defined by a stochastic nonlinear differential 
equation. Apart from direct numerical simulations the standard analytical tool as regards 
scaling properties is a perturbative renormalization grou p sc heme based on an expansion in 
powers of the nonlinear term in Eq. (jl.lj) or Eq. (jl.4j) jlll. IitI. llS^. This procedure yields 
renormalization group equations in an e-expansion about d = 2, e = d — 2, predicts a kinetic 
phase transition above d = 2 from a smooth Edwards- Wilkinson phase (£ = (2 — d)/2, z = 2) 
to a rough Burgers-KPZ phase with nontrivial exponents. In d = 1 the scheme yields 
(fortuitously) the exact exponents ( = 1/2 and z = 3/2. The limitation of the Langevin 



description is that it does not provide a simple physical picture of a growing interface and 
that the role of the noise can only be discussed and interpreted on a qualitative level. 

On the Fokker-Planck level the growing interface is determined by the deterministic evo- 
lution equation (J1.23j) driven by the Hamiltonian (|1.15|) . The formal structure of Eq. (J1.23j) 
is equivalent to a functional Schrodinger equation in Wick-rotated imaginary time and allows 
via the Martin-Siggia-Rose functional integral for a mapping of the growth problem onto 
a non-Hermitian quantum field theory as discussed in paper II, see also [4]. The quantum 
field formulation, in addition to also providing an alternative framework for perturbative 
dynamical renormalization group theory following for example the Callen-Symanzik scheme 
371 ] , permits two new lines of approach to the growth problem. Firstly, by a mapping of the 
Martin-Siggia-Rose path integral onto a directed polymer in a quenched random medium 



18, 



221 ] the nonequilibrium problem becomes equivalent to a disorder problem affording a 
different perspective on the growth problem and yielding new insight. The second line of 
approach which we adhere to in the present context is to discuss the non equilibrium prob- 
lem directly in terms of field theoretical constructs. The original stochastic fluctuations on 
the Langevin level are then interpreted as quantum fluctuations on the Fokker-Planck level, 
where the noise strength A in Eqs. (jl.23|) serves the role of an effective Planck constant. 

In the context of canonical quantization the quantum field theory or quantum many 
body theory for the interface is defined by Eq. (jl.l5|) with the canonical momentum p = 
A5/5u replaced by the momentum operator p = —iA5/5u in a u-diagonal basis obeying 
the canonical commutation relation [p(x),u(x)] = — iA5(x — x'). In the Edwards- Wilkinson 
case for A = we are dealing with a free field theory and the elementary excitations or 
(undressed) quasi-particles are the linear non-propagating diffusive modes with quadratic 
dispersion uj = uk 2 , yielding according to spectral properties the dynamic exponent z = 2, 
defining the Edwards- Wilkinson universality class. As discussed in paper II it is also an 
easy task to evaluate e.g., the slope correlations (jl.l4j) as a purely quantum many body 
calculation. In the nonlinear Burgers case for A / we obtain correction to the linear mode 
dispersion law. Moreover, the quasi-classical analysis for A — > in papers II and III also 
identifies a nonlinear soliton excitation with dispersion law E oc A^ -1 / 2 !! 3 / 2 . A detailed 
analysis of the non-Hermitian quantum field theory has, however, not yet been achieved and 
will be considered elsewhere. 

In the present paper we continue our investigation of the noisy Burgers equation for the 
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nonequilibrium growth of an interface. We make use of the weak noise canonical phase space 
approach developed in paper III and consider the following important issues: i) The detailed 
growth morphology based on the multi-soliton many body description, ii) the non-Hermitian 
properties of the superimposed linear mode spectrum and the phenomenon of dynamical 
pinning and mode transmutation, and iii) the correlations in the Edwards- Wilkinson case, 
the anomalous diffusion of growth modes and switching and pathways in the Burgers-KPZ 
case, correlations in the multi-soliton sector, and correlations and scaling in the two-soliton 
sector. With respect to i) we stress that one of the advantages of the quasi-classical weak 
noise phase space approach propounded in papers II and III is that it provides a many body 
description of a growing interface in terms of solitons and linear modes. The Landau quasi- 
particle picture emerging on the Fokker-Planck level was discussed heuristically in paper 
II. Here we analyze in more detail the time evolution of a growing interface in terms of its 
elementary excitations. Regarding ii) we note that superimposed on the nonlinear solitons 
are linear modes obtained by a linear stability analysis of Eqs. (|1.6J) and (jl.7|) about a 
soliton mode. An analysis of the linear mode spectrum was initiated in paper I and II. In 
the present paper we complete the analysis also for a multi-soliton state and demonstrate 
among other properties that the linear modes subject to the nonlinear soliton modes undergo 
a mode transmutation from diffusive non-propagating behavior in the absence of solitons 
to propagating behavior in the soliton case. Finally, with regard to iii) on the scaling 
properties of the slope correlations ()1.14j) we provided in paper II only a heuristic expression 
for the scaling function F based on a general spectral representation. Here we amend this 
situation and present an explicit expression for F within the two-soliton approximation, 
or brief accounts of the present work we refer to (34], 3^; moreover, the papers [3^ and 



4Cj present a numerical analysis of the soliton-bearing mean field equations and a tutorial 
review, respectively. 

The paper is organized in the following manner. In Sec.|H]we discuss the growing interface 
in terms of soliton modes, in Sec. IIHI we consider the superimposed linear mode spectrum 
and discuss the mode transmutation alluded to above, in Sec. II VI we address the statistical 
properties, and consider the anomalous diffusion of growth modes, switching and pathways, 
correlations in the multi-soliton sector, and in detail the correlations in the tractable two- 
soliton sector. In Sec. |V|we present a summary and a conclusion. 
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II. A GROWING INTERFACE 



A growing interface governed by the noisy Burgers equation is a simple prototype of 
an intrinsically open and driven nonequilibrium system. In the noiseless case for rj = the 
interface is damped and the slope field u evolves subject to a transient pattern formation 
consisting of propagating and merging right hand solitons connected by ramp solutions, 
with superimposed damped linear modes l[ 41]. The motion is deterministic and non- 
fluctuational. At long times the solitons die out on a time scale set by 1/uk 2 , where k is 
the diffusive mode wavenumber. In the noisy case for rj ^ the time evolution and pattern 
formation change. The noise balances the damping and drives after a transient period the 
slope field into a stationary morphology composed of amplitude-matched right and left hand 
solitons with superimposed linear modes. The noise strength A is an essential parameter 
changing the qualitative morphology of an interface; this is reflected mathematically in Eq. 
(jl.llj) for the transition probability which has an essential singularity for A — ► 0. 

The above behavior is illustrated in the linear case where an explicit solution of Eq. (jl.lj) 
for a wavenumber component u^, Uk(t) = J dx u(x, t) exp(— ikx), of the slope field driven by 
the noise wavenumber component 1]^ is given by 

Uh {t) = u\e-" kt - f dt'ike-^-^rikit'). (2.1) 



Here u>k = vk 2 is the diffusive mode dispersion law and u l k = Uk{t = 0) the initial slope value. 
We notice that generally l/cu/c sets the time scale. Initially the motion is deterministic and 
governed by the noiseless diffusion equation; at longer times for Lo^t ^> 1 the noise gradually 
picks up the motion as indicated by the kernel exp[— Uk(t — t')] in Eq. ()2.1|) and uu begins 
to fluctuate and is driven into a stationary noisy state. This behavior is in accordance with 
the phase space behavior discussed in paper III on the Fokker-Planck level. In Fig. ^ we 
have for a particular noise realization depicted the behavior of Uk- We emphasize that the 
general aspects of the noise-induced time evolution also holds in the noisy Burgers case here 
subject to a soliton-induced pattern formation. The transient regime is indicated by I, the 
long time stationary regime by II. 

As mentioned above the quantum mechanical interpretation allows a discussion of the 
growing interface in terms of a Landau quasi-particle picture. In the Edwards- Wilkinson 
or noninteracting case the relevant quasi-particle is the diffusive mode Uk with quadratic 
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dispersion law uj = vk 2 . In the Burgers case it is a general feature of the Landau quasi- 
particle picture that interactions usually give rise to a dressing effect of the free (bare) 
quasi-particle, e. g., the inducement of an effective mass. In the diffusive mode case this 
corresponds to a dressing of the damping constant v leaving the dynamical exponent z = 2 
in u = vk 2 unaltered. However, as shown in papers II and III even for weak A a new quasi- 
particle emerges, the nonlinear soliton excitation, with dispersion law to oc (X/u 1 ^ 2 )^^ 2 . 

From a heuristic point of view we can regard the soliton as a self-bound state of diffusive 
modes; in other words, the solitons condense or nucleate out of the diffusive mode field. We 
note, however, that while the formation of localized soliton modes with superimposed linear 
modes is a well-known feature of deterministic evolution equations, e.g, the sine-Gordon 
equation and the nonlinear Schrodinger equation 4^|, the underlying mechanism of the 
doubling of soliton modes here is the noise. In the approach in jjj the soliton mode was 
identified as a special solution of the classical field equations (jl.6|) and ()1.7|) obtained in 
the limit A — > from i) in |4| the underlying Heisenberg field equations pertaining to the 
quantum description and ii) in papers II and III from the classical field equations arising 
from a principle of least action in the WKB limit of the Fokker-Planck description. Below 
we turn to a discussion of the fluctuating interface in terms of the quantum/classical picture 
discussed above. 



A. The Edwards-Wilkinson case - equilibrium interface 



In the Edwards- Wilkinson case [23, 35] the slope of a fluctuating interface is governed by 
the driven conserved diffusion equation 

^ = isV 2 u + Vt 1j (2.2) 

which is readily solved both on the Langevin level in Eq. (|2.1j) and on the Fokker-Planck 
level. Since the damping term u'V 2 u in Eq. (j2.2j) can be derived from a thermodynamic free 
energy F = (1/2) J dx u 2 the driven diffusive equation describes an interface in equilibrium 
at a temperature T = A/2u. In accordance with the quantum field interpretation outlined 
above, this simple case, however, serves as an illustration of the quasi-particle representation. 
Consequently, we turn to the field equations (jl.6|) and (jl.7|) in the linear case for A = 0: 

^ =vV 2 u-V 2 p, (2.3) 
10 



% =-^VV (2.4) 



For a single wavenumber component Eq. ()2.2|) corresponds to a noise-driven overdamped 
oscillator with force constant u k = vk 2 and the associated canonical field equations ()2.3j) 
and ()2.4|) were solved and discussed in paper III. We find, supplementing the analysis in 
paper III, the solutions 

t4 sinh c^t + u\smhuj k (T - t) 
Uk{t) = sinh^T ' (2 ' 5) 

Pk® =^ Uk . Uk \ , (2.6) 
smh uj k l 

for an orbit from u\ to u k in time T, < t < T. The noise field Pk{t) is slaved to the 
motion of u k and determined by u k , u f k , and T. During the time evolution it evolves from 
p k — v[u{ — u k exp(— u k T)]/ smh u k T to p[ = u[u{ exp(u; fc T) — u),]/ sinhu; fe T. 

We note the correspondence between the physical interpretation on the Langevin level 
given by Eq. (|2.1j) and on the Fokker-Planck level characterized by Eqs. (J2.5)) and (J2.6|) . 
In the noiseless case for 77 = 0, corresponding to setting p k = 0, i.e., u\ = u k exp(—UkT), 
the slope field is damped according to Ukif) = u k exp(—ujkt) over a time scale l/u>k- In 
the presence of noise the growing noise field pk oc exp(ukt) eventually drives u k , i.e., Uk oc 
exp(uikt). Generally, Uk is a linear combination of a damped part exp(— u^t) and a growing 
part exp(djfct), analogous to the decomposition of the field in positive and negative frequency 
parts in quantum many body theory [43]. Here the components are decaying and growing 
according to the transient and stationary regimes I and II in Fig. Q respectively. 

The orbit (v,k,Pk) given by Eqs. ()2.5|) and ()2.6|) . representing the quasi-particles in the 
classical limit A — > 0, is confined to a submanifold in phase space delimited by four global 
conservation laws: Conservation of energy E = H, conservation of momentum II, conserva- 
tion of area, i.e., the integrated slope M or height offset; and here also conservation of the 
integrated noise field P given by Eqs. (|1.15|) . (jl.l7j) .and (jl.l8|) for A = 0, respectively. 

In wavenumber space we have 

E = E k = k 2 p. k (p k - 2uu k ), (2.7) 

J 2tt J Ait 

and the energy decomposes in contributions E k for each wavenumber mode k. Inserting Eqs. 
()2.5|) and ()2.6|) we obtain specifically 

uj 2 k \uj\ 2 + |4 1 2 - 2u{uj k cosh^T 

h k = T-T2 7f > \ 2 - H ) 

2 smh uj k l 
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and the energy only depends on the initial and final configurations u k , u\ and the time 
interval T. For fixed u k and u k and in the long time limit T —>■ oo the energy E k — > and 
the orbit migrates to the zero energy manifolds: pk = 0, the transient noiseless submanifold, 
and pk = 2uuk, the stationary noisy submanifold. The orbit thus asymptotically passes 
through the saddle point (uk,Pk) = (0,0) where the diverging waiting time ensures ergodic 
behavior. In Fig. El we have depicted the orbits in (u k ,Pk) phase space. 
In a similar manner, the momentum II decomposes according to 

„ f dk „ r dk , / > 

II = / — Il k = / — iku.kVk- (2.9) 

J ATI J ATI 

We note that II vanishes on the zero-energy manifolds pk = and pk = 2uuk', in the latter 
case since the integral in Eq. (|2.9J1 becomes a total derivative. For a finite time orbit 
insertion of Eqs. ()2.5|) and ()2.6|) explicitly yields 

* = (2.io) 

sinn uj ki- 
rn, terms of u k , u k , and T; for T — ► oo we have Hfc — > 0. Likewise for the integrated slope 
and noise field, P = 2uM, we have 

M = J dx u\x) = u k=0 . (2.11) 

Finally, the action associated with an orbit is obtained from Eq. Ijl.lOj) . Inserting the 
equation of motion Eq. ()2.3|) we have as an intermediate result S = (l/2iv) J dkSk, Sk = 
(1/2) / dtk 2 \p k \ 2 and using Eq. ()2.6j) the action 

s = u f dk l4~4 ex P(-^T)| 2 

J 2tt 1 - exp(-2o; fc T) ' 1 ' ' 

determined by the initial and final configurations u\ and u\ and the traversal time T. 
According to Eq. (jl.llj) we subsequently obtain the transition probability 

v f dk \u{ — w 1 fc exp(— u k T)\ 2 



p (u k ,ul,T) oc exp 



(2.13) 



A J 2tt 1 - exp(-2cu fe T) 

a well-known result |44i l45j ) . In the limit T — > oo the orbit migrates to transient-stationary 
zero-energy manifolds and we arrive at the stationary Gaussian distribution (jl. 21)1 . This 
behavior in phase space is consistent with the qualitative behavior shown in Fig. ^ 

Summarizing, in the linear Edwards- Wilkinson case the conserved noise-driven stochastic 
diffusion equation is in the weak noise limit equivalent to coupled field equations admitting 
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both damped and growing solutions for the slope field. The stochastic noise is replaced by 
a noise field canonically conjugate to the slope field. Both damped and growing solutions 
are required in order to describe the crossover from the transient regime to the stationary 
regime. The wavenumber k is a good quantum number and we can envisage the fluctuating 
interface as a gas of growing and damped diffusive modes according to the decomposition, 
see also paper II, 



In field theoretical terms u(x, t) is a free field and the elementary modes are noninteracting. 
A particular mode lies on the energy surface and is moreover specified by the conserved 
momentum I!*,. Furthermore, under time evolution the integrated slope field M and noise 
field P are also conserved. Finally, with the mode is associated an action yielding the 
transition probability P. 

The description based on the field equations (|2.Hj) and ()2.4j) and the associated symplectic 
structure is basically classical. Subject to canonical quantization the diffusive modes become 
bona fide elementary excitations and a Landau quasi-particle picture of the interface emerges. 
The original noise fluctuations are then interpreted as quantum fluctuations emerging from 
the underlying operator algebra. Finally, we note that subject to a Wick rotation t — > it the 
diffusive quasi-particles are transformed to dispersive particle-like quasi-particles with mass 



B. The Burgers-KPZ case - nonequilibrium interface 

In the case of a growing interface the situation is more complex. The general behavior 
depicted in Fig. ^ still holds in the sense that the interface evolves from a transient state 
to a stochastic stationary state. However, unlike the linear equilibrium case where the 
fluctuations in u are extended and diffusive, the fluctuations in the nonlinear nonequilibrium 
growth case include localized propagating modes in order to account for a growing height 
profile. This is also evident from e.g., the KPZ equation in Eq. (jl.4j) . where the damping 
and growth terms transform differently under time reversal. For t — > — t, h — > —h, and 
v — > —v the KPZ equation stays invariant. This is consistent with the fact that in the 
decomposition of the irreversible linear modes in growing and decaying components the 





A/2z/. 
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damping v enters in the combination ut, whereas the average nonlinear reversible growth 
term in the stationary state (dh/dt) st = (X/2)((Vh) 2 ) is invariant. Clearly, the growth 



term cannot be derived 
potential condition 



iom a thermodynamic free energy, the term moreover violates the 
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48, 



49| and drives the system away from thermal equilibrium 



into a stationary kinetic growing state. 

The issue on the Fokker-Planck level is again to establish a quasi-particle picture and to 
determine orbits in (u, p) phase space from an initial configuration u l at time t = to a final 
configuration u l at time t = T in order to, via the action associated with the orbit, evaluate 
the transition probability P(u { , u 1 , T). The orbit is in principle determined as an initial-final 
value problem, i.e., a boundary value problem in time, of the mean field equations (jl.6|) and 
(jl.7|) . Unlike the linear case where we can expand u on plane wave diffusive modes and thus 
achieve a complete analysis, the nonlinear and presumably nonintegrable character of Eqs. 
(II. 6|) and (jl.7j) precludes such an analysis. 

Two new features distinguish the field equations (jl.6|) and (jl.7|) from the linear case: i) 
the nonlinear coupling strength A setting together with v an intrinsic length scale vj\ and ii) 
the amplitude- dependent Galilean invariance (jl.Hj) . The new length scale allows for the pos- 
sibility of localized nonlinear excitations and the Galilean symmetry permits the generation 
of a class of propagating particle- like excitations from a static solution, i.e., an excitation at 
rest. The static excitations are the right and left hand solitons ()1.8|) characterized by the 
parity index fi = ±1. Boosting a static soliton to the velocity v, denoting the boundary 
values u + and w_, and using the Galilean symmetry in ()Oj) we obtain the fundamental 
soliton condition 

2v 

u+ + u_ = — -. (2.15) 
A 

In Fig. El we have depicted the fundamental "quarks" or solitons and the associated height 
profiles. 

The right hand soliton for = +1 moves on the noiseless transient sub manifold p = 
and is a solution of the damped noiseless Burgers equation, i.e., Eq. (II. 1J) for r] = 0. Within 
the canonical framework the right hand soliton does not contribute to the dynamics of the 
interface; according to Eqs. (|1.15|) and (jl.lfij) with p — it carries zero energy and zero 
momentum. The left hand soliton for fi = —1 is associated with the noisy or stationary 
submanifold p = 2uu and it follows from the field equations that it is a solution of the 
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growing noiseless Burgers equation for v — > —v. Note that due to the uneven boundary 
values u + and w_ the solitons are self-sustained dissipative structures driven by boundary 
currents as discussed in paper I. The left hand soliton is endowed with dynamical attributes 
and carries according to Eqs. (|1.15|) . and (jl.lUJl energy, momentum, and action: 

E = ^uX(u 3 + -u 3 _), (2.16) 

n =v{u 2 + -u 2 _), (2.17) 

S = -uX\u+ -m_| 3 T. (2.18) 
6 

Since u + < for a left hand soliton its energy is negative. Expressing II in the form 
II = — (2vi//\)(u+ — U-) using Eq. ()2.15|) it follows that II points in the direction of v. From 
II = mv we can also associate an amplitude-dependent mass m = (2u/X)\u + — U-\ with the 
soliton. Finally, the action for a left hand soliton orbit over time T is positive and Galilean 
invariant. 

In addition to the localized soliton modes the field equations also support linear modes 
superimposed on the soliton. These modes are obtained by a linear stability analysis of the 
field equations and will be discussed in Sec. IIHI . In the limit A — > the soliton modes vanish 
and the remaining fluctuations are the diffusive modes of the Edwards- Wilkinson model. 

The field equations (jl.fi)) and (| 1-7)1 are nonlinear and the soliton solutions do not con- 
stitute a complete set in the same manner as the plane wave decomposition ()2.14)) in the 
Edwards- Wilkinson case. We shall nevertheless as a working hypothesis assume that we 
can resolve a given initial interface slope profile u in terms of a gas of right hand and left 
hand solitons matched according to the soliton condition ([2.15)1 : i.e., with horizontal con- 
stant slope segments. Associated with the soliton representation of u there is also a soliton 
representation of the associated noise field p. From the form of the field equations it follows 
that a multi-soliton configuration is an approximate solution provided we can control the 
overlap contribution Xu± Vw w arising from the nonlinear term; here % and I denote the i-th 
and 1-th soliton. Since the soliton width is of order u/Xu the overlap only contributes in a 
region of order vjXu about the soliton center and is small in the inviscid limit v — > and 
for a dilute soliton gas. Otherwise, we assume that at least for small v we can absorb the 
correction term in a linear mode contribution. In summary, we represent a slope configura- 
tion u and the associated noise field configuration p in terms of a gas of right hand and left 
hand solitons matched according to Eq. ()2.15)) with superimposed linear modes. 
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This representation of a growing interface is in t 



43, 



re sp irit of a Landau quasi-particle picture 
5(1 . In a heuristic sense we assume that 



of an interacting quantum many body system 
the interface at a given instant of time is characterized by a gas of localized soliton modes 
and extended linear modes. Since the soliton is not associated with a particle but is a 
nonlinear solution of classical field equations the soliton number is not conserved; in other 
words, solitons are created and annihilated subject to collisions. 

Keeping only the solitons, explicit expressions for u and p are easily constructed in terms 
of the Galilee-boosted soliton modes ()1.8|) . Introducing the mean amplitude k p and velocity 
v p in terms of the boundary value u p+ i and u v of the p-th soliton, 

k p = -^(u p+ i - u p ), (2.19) 
v p = ~{u p+1 + Up), (2.20) 

we obtain for an n-soliton representation of the slope field u s and associated noise field p s 

2v n 

u s (x,t) — — fcptanh \k p \(x — v p t — x p ), (2.21) 
A p=i 
(2z/) 2 n 

p s (x,t) = — - — ^ fcptanh \k p \(x — v p t — x p ). (2.22) 

^ p=l,fc p <0 

The solitons are arranged from left to right according to the increasing index p, p = 1, 2, ■ ■ • n. 
The center of the p-th soliton is at x p and we have set U\ = u n+ i = 0. In Fig. HJwe have 
shown an n-soliton configuration. 

The soliton representation in Eq. (j2.21j) of the interface evolves in time according to 
the field equations. The motion corresponds to an orbit in (u,p) phase space lying on the 
manifold determined by the conservation of energy E, momentum II, and area M according 
to Eqs. 15j) . ()1.16j) . and (|1.17|) . We also observe from Eq. (jl.lSJ) that the integrated noise 
field P, since II is conserved, develops linearly in time according to 

P = XUt + const. (2.23) 

In the soliton representation the contributions to the energy, momentum, action, and area 
decompose. Noting that only left hand solitons for u p+ i < u p contribute dynamically we 
have applying Eqs. (|2~TH)l . (l2~T7l) . (I2~TH1) and (flTTI) 

E = -u\ Yl (4+1-4), (2.24) 

P=l,Wp+l<Up 
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n 



n = v 



E K 




(2.25) 



1 n 
S = -v\T V 




(2.26) 



■n 



p=i 



(2.27) 



Although, as follows from Eqs. (J2~24j) . (f2~25jl . and (|2~2fij) . the total energy, momentum, 
and action are additive quantities (extensive) the soliton gas still represents a very intricate 
many body problem. This is due to the soliton matching condition in Eq. (J2.15|) . i.e., the 
horizontal segments connecting the solitons, and the dynamical asymmetry between left 
hand and right hand solitons. The solitons in the representation (|2.21|) propagate with 
in general different velocities v p and are thus subject to collisions. Since we only have at 
our disposal single soliton solutions of the field equations patched together to represent a 
slope configuration at a particular time instant and not a general solution we have limited 
control over soliton-soliton scattering. Clearly, the expressions ()2.24l - 12.26)1 only hold in- 
between soliton collisions. In particular, the time T entering in the action ()2.26|) refers to 
times between collisions, i.e., T is typically of order \x p+ i —x p \/v p . The working assumption 
here is that in-between collisions the soliton (plus linear modes) representation is valid and 
that energy, momentum, and area are conserved during collisions. However, the number 
of solitons is not preserved, i.e., solitons are created and annihilated subject to collisions. 
Finally, we note that at long times the orbit from u 1 to w f migrates to the zero-energy 
manifold as conjectured in paper III. This implies that the finite energy solitons at long times 
are suppressed and that the system in this limit is described by diffusive modes yielding the 
stationary distribution (jl.21|) . For further illustration we have in Fig. depicted the slope 
field u, height field h, and noise field p for a 4-soliton configuration. The solitons are centered 
at X\, X2, x 3 , and X4 and propagate with velocities v\ = — (A/2)w 2 , ^2 = — (^/2)(m3 + M2), 
^3 — — (A/2)(m 4 + M3), and w 4 = —(X/2)u 2 , where u 2 , u 3 , and w 4 are the plateau values. The 
configuration carries energy E = (2/3)^A(w2 ~ "^3)5 momentum II = v{u\ — M3), action S = 
(1/6)^AT(|m 2 | 3 + |u4— M3| 3 +|«4| 3 ), and area M = U2(x2~Xi)+u 3 (x 3 —x 2 )+u 4: (x4 : —x 3 ) at time 
t = 0. By integration we note that the total noise field P evolves like P = Pq + \vt(u\ — u 3 ) 
in agreement with Eq. ()2.23j) . 

It is clear that the nonequilibrium growth of the interface is fundamentally related to the 
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existence of localized propagating soliton modes. Expressing the KPZ equation ()1.4|) in the 
form 

f = + \u 2 + V} (2.28) 

the linear damping term is associated with the linear modes, whereas the nonequilibrium 
growth term is driven by the solitons. In the Edwards- Wilkinson case the modes are extended 
and diffusive, u oc J2k fk(t) cos(kx + ip), where fk(t) is related to Eq. ()2.5j) . and the height 
field h = J x udx' for a particular k-mode behaves like h oc sm(kx + <£>); consequently, (h) = 
and growth is absent. On the contrary, in the Burgers case the localized soliton modes 
emerge and h(x) = J x u(x')dx' grows owing to the propagation of solitons across the system. 
Averaging Eq. ()2.28|) and setting v = we have (h) = (A / '2) {u 2 )t + const, which is consistent 
with the passage of a soliton with amplitude u and velocity \v\ = (X/2)u at a given point x. 
This growth behavior also follows from inspection of Fig. EJ 

The constant slope and noise field configuration u = uq and p = po have according to 
Eqs. (II. 15(1 and (jl.l6|) vanishing energy E and momentum II and thus form a continuum 
of zero-energy states; note that the energy is not bounded from below, i.e., the zero-energy 
states are not ground states. The right and left hand solitons in Eq. (II. 8|) lift the zero-energy 
degeneracy and connect a constant w_ configuration to a constant u + configuration. Unlike 
the solitons in the ip 4 theory or sine-Gordon theory J3, |^, connecting two degenerate 
ground states ±y? or degenerate ground states ip = irp, where p is an integer, respectively, 
with massive gapful excitations, the Burgers solitons are gapless modes forming a continuum 
of states with energy E oc (u 3 ^ — u^_) and momentum II oc (u 2 + — u 2 _). 

In the discussion of a growing interface in terms of its slope field we must introduce 
appropriate boundary conditions in order to describe the physical growth state. In the 
instantaneous configuration in Figs. E] and El we chose for convenience vanishing slope at 
the boundaries. However, owing to the soliton propagation this boundary condition cannot 
be maintained in the course of time as the solitons cross the boundaries of the system 
and it is more appropriate to assume periodic boundary conditions for the slope field, i.e., 
u(x) = u(x + L) at all times, where L is the size of the system. Note that periodic boundary 
conditions for the slope field in general does not imply periodic boundary conditions for 
the associated height field h, the integrated slope field. We have from h(x) = J x u(x')dx' , 
h(x) = h(x + L) + M, where M is the area under u, and only in the case of zero-area slope 
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configurations does h also satisfy periodic boundary conditions, corresponding to vanishing 
height offset at the boundaries. 

Whereas periodic boundary conditions for the slope field are consistent with the extended 
diffusive modes in the Edwards- Wilkinson case, i.e., the free fields, the elementary right and 
left hand Burgers solitons violate the boundary conditions since they connect unequal zero- 
energy states. In this sense we can regard the solitons as "quarks" in the present many body 
formulation. A proper elementary excitation or quasi-particle satisfying periodic boundary 
conditions is thus composed of two or more "quarks" as illustrated in Figs. E| and El 



1. The two-soliton configuration 



The simplest configuration satisfying periodic boundary condition is composed of two 
solitons of opposite parity, i.e., a noisy and a noiseless kink. The solitons have the common 
amplitude u, are centered at X\ and x% and propagate as a composite entity, according to 
Eq. (|2.15jl with velocity v = —Xu/2. Specifically, this pair-soliton mode has the form 



k x 



k K 



tanh — (x — vt — x\) — tanh — (x — vt — x-i, 



2 v ' 2 

The configuration w 2 together with the associated noise field p 2 (for u > 0) 



(2.29) 



P2OM) = vu 



and the height field 



h 2 {x,t) 



it 



k 

1 — tanh — (x — vt 
2 v 



cosh — vt — Xi) 



X2 



(2.30) 



k s °^ cosh !f(x — vt — x 2 ) 



+ const. 



(2.31) 



are depicted at the initial time t = in Fig. |U1 

According to Eqs. ()2.24l - 12.27)) the two-soliton configuration is endowed with the dy- 
namical attributes: 



E 

n 



2 \ 1 ,3 

--v\\u\\ 
i/sign(w)w 2 



S = -v\\u\ 3 T, 
6 

M = u(x2 — Xl). 



(2.32) 
(2.33) 
(2.34) 
(2.35) 
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The two-soliton configuration corresponds to an orbit in (u,p) phase space. Choosing as 
initial configuration at t — (u l ,p l ) = (u2(x,0),p2(x,0)) the final configuration at t — T 
is then given by (w f ,p f ) = (u2(x 1 T),p 2 (x 1 T)). In a finite system of size L with periodic 
boundary conditions this is moreover a periodic orbit with period L/v; in the thermodynamic 
limit L — > oo the period diverges. 

By inspection of Fig. H it follows that the two-soliton configuration propagates with a 
constant profile preserving the area M, i.e., the height offset 2u(x2 — x\). Subject to periodic 
boundary conditions in a system of size L the soliton pair reappears after a period L/v. This 
motion corresponds to a simple growth scenario where a layer of thickness Ah = u(x 2 — x\) 
is added to h per revolution of the pair. Subject to this particular soliton mode the interface 
thus grows with velocity (l/2)Au 2 (x2 — Xi)/L. This is consistent with the averaged form of 
Eq. ()2.28|) in the stationary state, (dh/dt) = (A/2)(w 2 ), noting the spatial weight {x2 — xi)/L 
of the soliton pair in the average (dh/dt) over the interface. We also remark that the local 
increase in h, Ah = u£, owing to the passage of a soliton pair of size i = \x 2 — x\\ in time 
At = £/\v\, where \v\ = Xu/2, yields Ah/ At = (X/2)u 2 , again is in accordance with the 
averaged KPZ equation. Finally, the integrated noise field P = J dxp decreases linearly 
with time as the soliton pair revolves like P = Pq — Ai>u\u\t in agreement with Eq. ()2.23|) . 

According to v = —Xu/2 the velocity of the soliton pair is proportional to the amplitude 
u. Expressing the energy E and momentum II in terms of v we have E = — (16/3)(z//A 2 )|v| 3 
and II = 4(i> / X 2 )v\v\ characterizing the nonlinear excitation. Moreover, eliminating the 
velocity we arrive at the dispersion law 



The soliton pair is thus a gap less quasi-particle mode with exponent z — 3/2. As discussed 
in paper II a general spectral representation for the slope correlations allows us to make 
contact with the scaling form in Eq. (jl.20j) and identify the mode exponent z with the 
dynamic exponent. 

Within the present description the statistical weight of the soliton is determined by Eq. 
(11.21)1 . In the inviscid limit for small v we obtain for a pair of size i and amplitude u the 
normalized stationary distribution 



E 



4 A 



nf/ 2 . 



(2.36) 



3^/2 



P st (u,£) = n- 1 (L)ex P --uH , 



(2.37) 
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with normalization factor 



(£) = £ dl J du £P st (u, £). (2.39) 



Qst(L) = 2(ttA/z/) 1/2 L 1/2 . (2.38) 

The distribution is parameterized by the amplitude u and the size £. The normalization 
factor or partition function f2 st varies as L 1 / 2 and the distribution thus vanishes in the 
thermodynamic limit L — > oo, characteristic of a localized excitation. The mean size of a 
soliton pair is given by 

Inserting P st we obtain (£) = (1/3)L, i.e., the mean size of the pair scales with the system 
size. This behavior is characteristic of a spatially extended or loosely bound excitation and 
we can envisage the soliton pair as a "string excitation" connecting right and left hand 
solitons (the fundamental "quarks"). 

In a similar manner we can determine the transition probability associated with a soliton 
pair using Eq. Inserting S from Eq. (J2.34j) and normalizing we obtain 

v\ 



P sol (u,T) = n^(T)exp 



6A 1 1 



^2.40) 



where the normalization factor or dynamic partition function is given by 



3 





■6A 1" 


G) 


.7\T. 



1/3 



(2A1) 



Here the Gamma function T(z) = J Q exp(—t)t z 1 dt arises from the normalization of P so \, 
r(l/3) = 2.68174. 

Before turning to the linear fluctuation spectrum in the next section we wish to briefly 
review a recent numerical study of the field equations (jl.fjjl and (jl.7|) j^. The field equa- 
tions are of the diffusive-advective type with the characteristic feature that the equation 
for p admits exponentially growing solutions due to the negative diffusion coefficient thus 
rendering direct forward integration in time numerically unfeasible. In order to resolve this 
instability problem we developed a "time loop" integration procedure based on integrating 
the equation for u forward in time followed by an integration backward in time of the equa- 
tion for p. This numerical scheme thus requires an assignment of both initial and final (u,p) 
configurations and therefore mainly served as a tool to check whether a certain assignment 
actually constitutes a solution. 
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We investigated numerically three propagating soliton configurations: i) a propagating 
soliton pair, ii) two symmetrical solitons colliding with a static soliton, and hi) the collison 
of two symmetrical soliton pairs. Referring to 39| for details we summarize our findings 
below. We found in case i) that the pair-soliton in the inviscid limit for small v is a long lived 
excitation thus justifying the quasi-particle interpretation above. In case ii) we considered 
symmetrical solitons propagating towards and colliding with a static soliton at the center 
passing through the static soliton and reemerging with no phase shift or delay; this specific 
mode corresponds to filling in a dip and subsequently nucleating a tip at the same point in 
the height field. In case iii) we finally considered two symmetrical soliton pairs colliding and 
reappearing subject to a phase shift or delay where the incoming trailing solitons become 
the leading outgoing solitons; this mode correponds to filling in a trough and subsequent 
nucleation of a plateau in the height profile. It was characteristic of the soliton collisions in 
case ii) and iii) that the conservation of E, II, and M was satisfied during collision, a feature 
that seems to stabilize the integration. 



III. FLUCTUATIONS - MODE TRANSMUTATION 

The soliton spectrum discussed in the previous section is a fundamental signature of the 
nonlinear character of the Burgers equation and at the same time essential in accounting 
for the growth aspects of an interface. In the present section we address the fluctuation 
spectrum or linear mode spectrum superimposed on the soliton gas. In the linear Edwards- 
Wilkinson case discussed in Sec. |Il] the fluctuations exhaust the mode spectrum and have 
a diffusive character. In the nonlinear Burgers-KPZ case the fluctuations are superimposed 
on the soliton modes and become propagating. The fluctuation spectrum was discussed 
incompletely in the noiseless Burgers case in paper I and in the noisy case in paper II. Here 
we present a detailed analysis. A brief account of the work is given in [3] ; here we complete 
the analysis. 
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A. The noiseless Burgers case 



The noiseless Burgers equation is inferred from Eq. (jl.lj) by setting 77 = and also follows 
from the field equation (Jl.fcij) on the p = submanifold, 

^ - AuvJ u = uV 2 u. (3.1) 

This equation exhibits a transient pattern formation composed of right hand solitons con- 
nected by ramps with superimposed linear modes. A single right hand static soliton mode 
is given by Eq. ()1.8|) for /i = 1, i.e., u s = uta.nhk s x, k s = Xujlv. The soliton has amplitude 
u and width k~ l . A spectrum of moving solitons is then generated by the Galilean boost 
(jl.3|) : x — > x — Xuot, u — > u — Uq. 

In order to analyze the superimposed linear fluctuations we expand u about the soliton 
mode u s , u = u s + 5u. To linear order in 5u we obtain the equation of motion 

- A^vJ 5u = uV 2 5u + X(Vu s )5u. (3.2) 

In the asymptotic limit for large \x\ this equation is readily analyzed. Noting that 
u s — > Msign(x), Vu s — > 0, and searching for plane wave solutions of the form 5u oc 
exp(— Ekt) exp(ikx) we identify a spectrum of complex eigenvalues, 

Ek = vk 2 — iAu/csign(x), (3.3) 

showing the non-Hermitian character of Eq. (|3.2j) . Introducing the phase velocity v = Xu 
the imaginary part of the eigenvalue combines with the plane wave part and yields the 
propagating damped wave form 

5U OC e -^ 2 t e ik(x+vtsign(x))_ (3^) 

The soliton mode thus gives rise to a mode transmutation in the sense that the diffusive 
mode in the Edwards- Wilkinson case 5u cx exp(— uk 2 t) exp(ikx) is transmuted to a damped 
propagating mode in the Burgers case with a phase velocity v depending on the soliton 
amplitude u. The mode transmutation, of course, also follows from the Galilean invariance 
in Eq. (jl.3j) since a shift of the slope field to the soliton amplitude u corresponds to a 
transformation to the moving frame x — ► x — Xut. For large positive x the mode propagates 
to the left, for large negative x the propagation is to the right, i.e., the mode propagates 
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towards the soliton center which thus acts like a sink. The phenomena of mode transmutation 
has also been noted by Schiitz [52J in the case of the asymmetric exclusion model, a lattice 
version of the noisy Burgers equation, in the context of analyzing the shocks, corresponding 
to the solitons in the present context. 

As noted in paper I the analysis of Eq. (|3.2j) can be extended to the whole axis by 
introducing a nonuniform gauge function g determined by the soliton profile in Eq. (jl.8|) . 



g(x) = k s tanhfc s x , k s = Xu/2u. 

By means of g we can express Eq. (13. 2J) in the Schrodinger-like form 

d5u 



(3.5) 



dt 



D(g)5u, 



(3.6) 



where the operator D(g) in the quantum mechanical analogue is given by the Hamiltonian 

2 



D{g) = -u{V + g{x)Y + vkl 



1 - 



cosh k*x 



(3.7) 



This equation of motion describes in imaginary Wick-rotated time a particle moving in the 
potential —2/ cosh 2 k s x subject to the imaginary gauge field ig. Absorbing the gauge field 
by means of the gauge transformation 



U (g) = exp 



g(x) dx 



cosh 1 k s x, 



(3.8) 



and using the identity (V + gf = U ~(g)V 2 U (g)' 1 , the non-Hermitian equation of motion 
(J3.6|) takes the Hermitian form 

d5u 



dt 



D(0)6u, 



(3.9) 



where 



Su = U(g)- L 5u. 



(3.10) 



We also observe that the gauge transformation in Eq. ()3.8j) has the same form as the 
Cole-Hopf transformation |53|, |54| applied to the static soliton solution u si see also papers 
I-III. 

The presence of the gauge function changes the spatial behavior of the eigenmodes and 
thus their normalizability. The imposition of physical constraints such as spatially localized 
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modes or asymptotic plane wave modes obeying periodic boundary conditions consequently 
give rise to a complex eigenvalue spectrum. Since we can "gauge" the mode problem in 
Eq. (|3.6|) to the exactly solvable Schrodinger problem defined by Eq. (|3.9|) we are able to 
complete the analysis. 

Setting 5u oc exp(— fit), where fl is the frequency eigenvalue, the spectrum of the ensuing 
eigenvalue problem, 

D(0)6u = Sl6u, (3.11) 



associated with -D(O) can be analyzed analytically 55, 56]. It is composed of a localized zero 



frequency mode and a band of phase-shifted extended scattering modes with eigenvalue 

tt k = u(k 2 + k 2 ). (3.12) 

The eigenmodes have the form 

5u oc — — — , Q = 0, (3.13) 

cosh k s x 

5u oc exp(ikx)sk(x), flk = v{k 2 + k 2 ), (3-14) 

, . k ~\~ ike tcLiili kq$ , . 

s k {x) = ' ■ (3.15) 

The x-dependent s-matrix s^(x) = |sfc(x)| exp[^ fc (x)] gives rise to a spatial modulation 
|sfe(a?)| = [{k 2 + k 2 tanh 2 k s x)/(k 2 + k 2 )] 1 ^ 2 of the plane wave near the soliton center over a 
range k~ x and a phase shift 5k(x) = tan _1 ((/c s tanh k s x)/k) + ta.n~ 1 (k s /k). For x — ► — oo, 
Sk{%) — > 1; for x — > oo, Sk(x) — >■ (k + ik s )/{k — ik s ) = exp(i5k), and Sk(x) becomes the 
usual s-matrix. We note that the bound state solution and its zero-frequency eigenvalue is 
contained in the scattering solution as a pole in the s-matrix for k — > ik s . 

Inserting the gauge transformation in Eq. ()3.8|) we obtain for the zero-mode 5u oc 
cosh~ 2 (k s x) oc Vw s which thus corresponds to the translation or Goldstone mode associ- 
ated with the position of the soliton. For the extended states the gauge transformation 
U provides an envelope of range kj 1 , i.e., 5u oc exp(ikx)sk(x) cosh~ 1 (k s x). The complete 
solution of the mode equation (|3.2|) thus takes the form 

5u = /I +J2B k e- Ukt -^—s k (x), (3.16) 
cosh k s x ~^ cosh k s x 

expressing the fluctuations of the slope field about the static right hand soliton u s ; A and 
B k = B*_ k are expansion coefficients. The first term in Eq. ()3.16|) is the time independent 

25 



translation mode. The second term corresponds to a band of damped localized states with 
a gap ukg in the spectrum Q k = v(k 2 + k 2 s ). Moreover, the scattering modes are transparent 
and phase-shifted by 5 k = 2 tan -1 (k s /k), implying according to Levinson's theorem that the 
band is depleted by precisely one state corresponding to the translation mode. 

In order to make contact with the asymptotic analysis yielding the spectrum in Eq. (|3.3j) 
we observe that the fluctuations given by Eq. ()3.16|) do not exhaust the spectrum. Since the 
gauge factor U = cosh -1 k s x provides a fall-off envelope we can extend the set of solutions 
by an analytical continuation in the wavenumber k in the same manner as the translation 
mode is retrieved by setting k = ik s . Thus shifting k — > k + in, where \k\ < k s , we obtain 
by insertion the complex spectrum 

E k ,K = v{k 2 + k 2 s - k 2 ) + 2iukn, \k\ < k s , (3.17) 

and associated fluctuation modes 

Su = . + £ 5 Ke -(^-- 2 )yM-w)_ Sk+lK (x), (3.18) 
cosh k s x cosh k s x 

where B k ^ K = B*_ k K since 5u is real. 

The expressions f|3.17|) and (J3.18J) provide the complete analytically continued solution 

of the fluctuation spectrum about the static soliton compatible with the imposed physical 

boundary conditions. For k = ±k s we recover the spectrum (J3.3j) of right hand and left hand 

extended gapless modes propagating towards the soliton center with velocity v = 2uk s = Xu. 

For \k\ < k s we obtain a band of gapful modes with localized envelopes propagating with 

velocity 2z/k towards the soliton center. Finally, for k = the envelope is symmetric, the 

spectrum E^o is real, and the mode has no propagating component; for k = and k = k s 

we retrieve the time independent translation mode. With the exception of the extended 

mode for n = ±k s , the envelope modes for \k\ < k s are dynamically pinned to the soliton. 

This phenomenon of localization or dynamical pinning of the modes is associated with the 

complex spectrum in Eq. ()3.17|) resulting from the non-Hermitian character of the eigenvalue 

problem. In all cases the modes are damped with a damping constant given by v(k 2 -\-k 2 — Hi 2 ). 

We mention that a non-Hermitian eigenvalue spectrum is also encountered in the context 

of flux pinning and the transverse Meissner effect in high T c superconductors (otI . Io^ . Here 

the imaginary gauge field is uniform and is given by the applied transverse magnetic field; 

in the present case the gauge field is spatially varying and given by the nonlinear soliton 

excitations. 
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In Fig. |7| we have depicted the spectrum in a plot of the imaginary part of Ek, K versus its 
real part. In Fig. |H]we have shown the associated characteristic fluctuation mode patterns. 
Specifically, in order to obtain a real extended mode propagating towards the soliton center 
with velocity v = 2uk s Xu from the left and with velocity — v from the right with continuous 
derivative at x = 0, thus extending the asymptotic solution (|3.4|) to the whole axis we form 
an appropriate linear combination from Eq. (j3.18j) with equal weights according to the 
assignments ±k and k = ±k s . Ignoring here the modulation factor Sk+iK,(x), which is easily 
incorporated, we obtain 

_ uk 2 t e~ ksX cosk(x — vt) + e ksX cos k(x + vt) , , 

Oil OC 6 ■ . lo.-LJ) 

cos k s x 

This mode is depicted in Fig. El 



B. The noisy Burgers case 

In order to discuss the fluctuation spectrum in the noisy case we must address the coupled 
field equations (jl.6[) and (jl.7|) and expand the slope field u and noise field p about the soliton 
configurations, in the single soliton case given by Eqs. (jl.8|) and f)1.9j) and in the multi-soliton 
case by Eqs. (j2.21|) and ()2.22|) . As also resulting from the analysis in paper II it is convenient 
to make use of a symmetrical formulation and introduce the auxiliary noise field if by means 
of the shift 

p = u(u-(p). (3.20) 
The field equations then assume the symmetrical form 

^ - AnV^J u = uV 2 (p, (3.21) 

^ - XuVj = zA7V (3.22) 

The single soliton solution u 1 ^ is given by Eq. (jl.SJ) and the associated noise solution by 

(fg = fiu s + const. (3.23) 

Expanding about a general multi-soliton configuration (u a ,(p a ), where u s is given by Eq. 
(j2.21|) and cp a by (modulus a constant) 

2v n 

<p s (x, t) = — ^2 \k p \ tanh \k p \(x — v p t — x p ), (3.24) 
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by setting u = u s + 5u and <p = <p s + 5<p, we obtain the coupled linear equations of motion 

f - Aw s V^J 5u = uV 2 5(p + X(Vu s )8u, (3.25) 

- Xu s V^j 8<p = vV 2 5u + \(V<p s )8u, (3.26) 

determining the fluctuation spectrum of superimposed linear modes. 

The analysis proceeds as in the noiseless case. Referring to Eqs. (|2.21|) and ()3.24|) we note 
that in the inter-soliton matching regions of constant slope and noise fields Vu s = Vy2 s = 0. 
The Eqs. (|3.25|) and (|3.26|) then decouple as in the linear Edwards- Wilkinson case and setting 
u s = u and looking for solutions of the plane wave form 5u,5(p oc exp(— E^t) exp(ikx) we 
obtain 5u ± 5(p oc exp(— E^t) exp(ikx), i.e., Su oc [A exp(— E^t) + B exp(— E^t)) exp(ikx), 
where the non-Hermitian complex eigenvalue spectrum similar to the noiseless case is given 
by 

Ej*r = ±vk 2 - ivk, v = Xu. (3.27) 

The 5u mode (and likewise the 5(p mode) thus corresponds to a propagating wave with both 
a growing and decaying component, 

5u oc (Ae' vkH + Be vkH )e ik{x+vt) . (3.28) 

These aspects are consistent with the general phase space behavior depicted in Fig. |2l 
whereas the propagating aspect as in the noiseless case is the manifestation of a mode 
transmutation from diffusive modes in the Edwards- Wilkinson case to propagating modes 
in the Burgers case. As indicated in Fig. H] the linear mode propagates to the left for u > 
and to the right for u < 0; for u = the propagation velocity vanishes and we retrieve 
the diffusive modes in the Edwards- Wilkinson case. We note in particular that for a static 
right hand soliton with boundary values u± = ±it, equivalent to the noiseless case discussed 
above, the mode propagates towards the soliton center which acts like a sink; for a static 
noise-induced left hand soliton with boundary values u± = =fm the situation is reversed and 
the modes propagate away from the soliton which in this case plays the role of a source. 

In the soliton regions the slope and noise fields vary over a scale k~ l and we must address 
the equations (j3.25J) and ()3.26|) . Introducing the auxiliary variables 

5X ± = 8u±8(p, (3.29) 
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and the general gauge function g s defined by the slope profile 



9s(x,t) = ^-u s (x,t), (3.30) 



the equations ()3.25)) and (|3.26j) take the form 

d5X ± , „, , , , A 



<)f ±D(±g s )5X ± - -(y Ua ± V<p s )5X*, (3 31) 

where D(±g s ) is the "gauged" Schrodinger operator 

D(±g s ) = -u(V ± g s ) 2 + ^u 2 s - ±Vip„ (3.32) 

describing the motion of a particle in the soliton-defined potential (X 2 /4u)u 2 — (A/2)Vy? s 
subject to the gauge field g s . In the regions of constant slope and noise fields we have 
Vu s = Vy? s = 0, u s = u, g s = Xu/2u, D(±g s ) -> -u(V ± Xu/2vf + (X 2 /Av)u 2 , and 
searching for solutions of the form SX ± oc exp(— E^t) exp(ikx) we recover the spectrum in 
Eq. ()3.27|) and since 5u = 5X + + 5X~ the mode in Eq. ([3.28)1 . In the soliton regions we 
have V(p^ = /^Vm^, where \i = ±1 for the right and left hand solitons, respectively, and one 
of the equations ()3.32|) decouple driving the other equation parametrically. 

The analysis proceeds in a manner analogous to the noiseless case. Introducing the Cole- 
Hopf transformation 

U(x,t) = exp ^— J g s (x',t) dx'^j , (3.33) 

and using the identity (V + g s ) 2 = t/V 2 f/ _1 we arrive at the coupled Hermitian equations 

- = ±U ±1 D(0)U^ 1 5X ± - ^{Vu s ± Vcp s )6X*, (3.34) 

which are readily analyzed in terms of the spectrum of D(0) summarized in Eqs 1)3.11)) to 
1)3. 15)1 . The exponent or generator in the gauge transformation in Eq. ()3.33)) samples the 
area under the slope profile u s up to the point x. For x — > 00, U — > exp(— XM/2u), where 
M given by Eq. ()1.17j) is the conserved total area. In terms of the height field h, u — V/i, 
M = h(+L) — h(—L) for a finite system and thus equal to the height offset across a system 
of size L, i.e., a conserved quantity under growth. Inserting the soliton profile u s in Eq. 
()2.21j) the transformation U factorizes in contributions from the individual local solitons, 
i.e., 

n 

U(x, t) = [] U p (x, t) sigD ^ , U p (x, t) = cosh" 1 k p (x - v p t - x p ). (3.35) 
p=i 
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1. The single soliton case 



Since the formulation of the linear mode problem in terms of Eqs. ()3.25|) and ()3.26|) 
deriving from the field equations (J3.21j) and (J3.22j) and yielding Eqs. (J3.34j) is entirely 
Galilean invariant, we can in analyzing the p-th single soliton segment of the multi-soliton 
configuration in Eq. (j2.21|) without loss of generality boost the soliton to a rest frame 
with zero velocity. Thus shifting the slope field of the p-th soliton by the amount (u p+ i + 
Up) corresponding to the propagation velocity v p = — (\/2)((u p+ i + u p ), and assuming for 
convenience that x p = the soliton profile is given by in Eq. (jl.8j) . Hence we obtain 
Vy£ = /iV< and D(0) given by Eq. (EP2J) . Noting that U = Ug and V< = fiuk s U^ we 
find for the fluctuations 

5X* = U^5X ± , U p = cosh" 1 k s x, (3.36) 

the Hermitian mode equations 

_ = ±D(0)6X ± - vk 2 s {ii ± 1)5X T , (3.37) 

which decouple and are readily analyzed by expanding 8X ± on the eigenstates of -D(O). We 
obtain 

6X» = 2 <A^t±Bl_ + {A (,) e -,u k t + B M e »n k t f hx Sk(x) ^ (338) 
cosh k s x " coshA; s x 

5X-» = + uB^^e^j^sJx) coshfc s x, (3.39) 

vkj 

describing the fluctuations of the slope and noise fields 5u = (5X + + 5X~)/2 and cp = 
(SX + — SX~)/2 or Sp = v5X~ about the static right hand (/i = +1) and left hand (/i = — 1) 
solitons. Ap , Bq 1 ^, A k ^\ and B% are integration constants fixed by the initial conditions. 

The first terms in Eqs. (I3.38j) and 1)3.39)1 are associated with the soliton translation 
modes 5X^ M oc Vw^ oc cosh -2 k s x which propagate with constant momentum X^m oc 
Aq] we note that the soliton position and soliton momentum are canonically conjugate 
variables. The envelope modulated plane wave terms in Eqs. ()3.38j) and ()3.39)) represent the 
fluctuations about the soliton. Since the s-matrix Sk(x) —>■ (k+ik s )/(k—ik s ) = exp(i5k), = 
2 tan _1 (fc s /fc) for x — > oo the soliton induced potentials are transparent and the fluctuations 
pass through the soliton only subject to the phase shift 5k and a spatial modulation. We also 
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note that Levinson's theorem implies that the band is depleted by one mode corresponding 
to the translation mode 5^|. 

Confining the fluctuations to the noiseless transient submanifold p = we have Sp = 0, 
i.e., 5u = 5(p, and we obtain for the right hand soliton (/x = +1) 5X~ = implying = 

k 



and .£>[ = and thus 



5X+ = 26u = — % + Ak 6 ± Sk{X) , (3.40) 

cosh k s x cosh k s x 

in accordance with Eq. (|3.16|) in the noiseless Burgers case, i.e., a translation mode and a 
band of damped localized pinned modes. 

Likewise, on the noisy stationary submanifold p = 2vu we require 5p = 2u5u, i.e., 5u = 
Sip, and we obtain for the noise induced left hand soliton (/i = —1) 5X + = entailing 
= and B% = and we have 



SX- = 26u = 7 2 - + ~ k , . (3.41) 



B^ 1] A { ~ l) e nkt e ikx s k (x) 
cosh 2 k s x cosh k s x 

composed of a translation mode and localized modes. However, unlike the fluctuations about 

the noiseless right hand soliton which are damped, the fluctuations associated with the noisy 
left hand soliton are growing in time. This behavior is consistent with the phase space plot 
in Fig. Eland the Edwards- Wilkinson case discussed in Sec. ILT1 

In general there are also fluctuations perpendicular to the submanifolds and we are led 
to consider the coupled equations ()3.38|) and ()3.39|) . The fluctuations are modulated by 
the gauge factors cosh k s x and cosh -1 k s x. As in the noiseless case the spatial modulation 
of the plane wave form allows us to extend the spectrum by an analytical continuation 
in the wavenumber k and in this manner match the spectrum to the inter-soliton regions. 
In fact, noting that 5X ± according to ()3.34|) decouples for Vw s , V<p s — > in the inter- 
soliton regions, we obtain by setting k —>■ k ± ik s the shift Q k = Z7 (^ 2 +s) ~~ * u (k 2 ^ 2i/cfc s ) 
and exp(ikx) cosh 1 * 1 k s x const, and we achieve a matching to the extended propagating 
modes. The gauge transformation in Eq. (J3.35|) thus permits a complete analysis of the 
linear fluctuation spectrum about a multi-soliton configuration. 



2. The two-soliton case 



In order to illustrate how the fluctuation spectrum is established across the soliton con- 
figuration and how the matching is implemented we consider the fluctuations about the 
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two-soliton configuration in Eq. ()2.29j) with associated noise field 



<p 2 (x,t) 



k k 
tanh ~^{ x — vt — x\) + tanh -^-(x — vt — x 2 



The Hermitian linear mode equations are given by 
Cole-Hopf transformation 

A 



with tt s 



U 2 (x, t) = exp 



A f x 

— / u 2 (x',t) dx 
2v J 



(3.42) 

u 2 , <f s = V?2 and the 

(3.43) 



The soliton pair propagates with velocity v = —Am/2 so in order to render the gauge trans- 
formation time independent and thus facilitate the analysis we boost the configuration to a 
rest frame, u 2 — > u 2 — u/2. Inserting u 2 we then obtain the specific gauge transformation 



Uo(x) 



cosh • 



[x - x 2 ) 



(3.44) 



cosh y( x — x i) 

a special case of ()3.35|) . 

In order to simplify the discussion of the mode equations ()3.34|) with -D(O) given by ()3.32|) 
we introduce the notation u± = {u/2) tanh k s /2(x for the individual 

solitons contributing to u 2 and tp 2 , i.e., u 2 = u 1 —u 2 —u/2 and (p 2 = u 1 +u 2 . A simple estimate 
of the potential (\ 2 /Au)u 2 2 - (A/2)V^ 2 in (j3~^ then yields u(k s /2) 2 [l - V+- VL], where 
we have also introduced the notation V± = 2/ cosh 2 (k s /2)(x — x±), for the soliton-induced 
potentials. Moreover, (A/2)(Vw2 ± Vy2 2 ) = ±v(k s /2)V±, and we arrive at the two-soliton 
mode equations 

d5X ± 



dt 



±U± 1 D(0)U? 1 5X ± =f v(k s /2) 2 V±5X Z 



(3.45) 



with -D(O) given by 



D(0) = -vV 2 + v{k s /2) 2 {\ - V + 



V- 



(3.46) 



In the regions of constant slope field V± ~ and 5X decouple. To the right of the soliton 
pair for x 3> x\,x 2 we have U 2 oc exp(k s (xi — x 2 )/2) exp(k s x/2) and we obtain the envelope 
solutions 5X ± oc exp[±u(k 2 + (k s /2) 2 )t]U^ 1 exp(ikx) which are matched to the asymptotic 
plane wave solution by setting k — > k±ik s /2, yielding 5X ± oc exp(=pz/fc 2 t) exp(ik(x — uk s t)) , 
i.e., a mode propagating to the right with velocity vk s = Xu/2 = v. To the left for x <C xi, x 2 , 
we have U 2 oc exp(k s (x 2 — x\)/2) exp(k s x/2) and correspondingly the envelope solutions 
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8X ± oc exp[±(i/(/c 2 + (k s /2) 2 )t]U 2 1 exp(ikx). Matching these solutions to the plane wave 
solutions by the analytic continuation k — > k±ik s /2, we obtain the same result as above. We 
note that the change in U 2 across the pair is given by exp(— k s (x 2 — xi)) = exp(— (\/2b / )u(x 2 — 
xi)) ~ exp(— (X/2u)M2), where M 2 is the area enclosed by the soliton pair. In the region 
between the solitons for ii < a; < x 2 , U 2 ~ exp(k s (xi + x 2 )/2) exp(— k s x/2) or for k — ► 
k =1= ik s /2 the modes 5X ± oc exp(^fuk 2 t) exp(ik(x + uk s t)), corresponding to propagation 
to the left with velocity — uk s = —Am/2 = —v. In the soliton region near x + we have 
V- ~ 0. Ignoring the translation mode and phase shift effects we thus have from (J3.45)) 
5X~ oc exp(z/(A; 2 + {k s ) 2 )t)U 2 l exp(ikx). By insertion in Eq. (|3.45|) we note that V^U 2 = 
V 2 exp(—k s x/2) ~ and that 5X + thus is decoupled from 5X~ , yielding the solution 5X + oc 
exp(— vik 2 + (k s /2) 2 )t)U 2 exp(ikx). A similar analysis applies in the soliton region near x 2 . 

For the plane wave components alone we then obtain, interpolating to the whole axis and 
incorporating the s-matrices according to Eq. (|3.15jl . 

5X + ~ e -Kfc a +(fc./2) a )t Cosh(fc a /2)(x - x 2 ) u xl2eikx _ _ ^ (3 ^ 4?) 

cosh(A; s /2)(x — x\) 

5X- ~ ^^fflt ^^j^ - ^ -Wa^fr - - x a ). (3.48) 

cosh(A; s /2)(x — x 2 ) 

We thus we pick up the phase shift 5k = 2tan _1 (/c s //c) across each soliton, the plane 
wave components are moreover spatially modulated by the gauge transformation U 2 = 
exp[— (A/2z/) J u 2 dx] sampling the area under u 2 . 

Finally, we boost the mode to the velocity v = —vk s = —\u/2v, shift the wavenumber 
k — ¥ k±ik s /2 for 5X ± and obtain the modulated plane wave associated with the propagating 
two-soliton mode with vanishing boundary conditions 



rv , „ fc 2.cosh(£; s /2)(£ - vt - x 2 ) ikx , , 

<5X + ~ e — ; 7 ^e lfcx s fc+ifcs/2 (a; - irf - Si)s fc+ifc . /2 (x -vt- x 2 ), (3.49) 

cosh(/c s /2 i )(x — vt — Xi) 



The interpretation of this result is straightforward. In the regions away form the pair-soliton 
we obtain plane wave modes with both a growing and decaying time behavior. Across the 
propagating two-soliton configuration the plane wave amplitude and form is modified by the 
gauge factor and the s-matrix. 



33 



3. The multi-soliton case 



In the multi-soliton case the slope configuration u s and the associated noise field (p s are 
given by Eqs. 1)2.21(1 and ()3.24() . The linear mode problem is defined by Eqs. ()3.34() with 
the gauge transformation U given by Eq. (|3.33|) . As discussed previously the extended 
plane wave modes in the inter-soliton regions connecting the solitons are subject to a mode 
transmutation to propagating waves with spectrum given by Eqs. (13.27(1 . In the soliton 
regions the analysis follows from a generalization of the single and two-soliton cases discussed 
above. 

A complete analysis is achieved by first noting that U also can be expressed in the form 



U(x, t) = exp 



~ fdt' (^u 2 s (x,t') + vVip s (x,t') 



(3.51) 



derived by differentiating U in Eq. (|3.33() . using the equation of motion in Eq. (|3.21() for 
the multi-soliton profile and solving the ensuing differential equation. By insertion of Eq. 
()3.51j) in Eqs. (|3.34() we obtain the linear equations of motion 

^ (u^SX*) = ±(z/V 2 + \V<p s )U Tl 5X ± + ^(Vu s ± V Vs )U Tl 5X T , (3.52) 

which are readily discussed. In the flat regions Vw s = V(p s = , 5X + and 5X~ decouple 
and we have 

5X ± oc e^ vkH e kx \J ±x [x, t), (3.53) 

describing a plane wave mode modulated across the soliton regions by the Cole-Hopf trans- 
formation U(x,t) with the explicit form given by Eq. ()3.35() . Alternatively, we obtain a 
mode transmutation to a propagating plane wave between solitons by the analytical con- 
tinuation k — ► k ± ik p thus absorbing the spatial modulation in U (x, t) and corroborating 
the previous discussion. We note that the form in Eq. (|3.53() is in accordance with Eqs. 
(13.49(1 and (|3.50J1 in the two-soliton case. Across the soliton regions SX + and 5X~ couple 
according to Eqs. ()3.52() . The analysis in the single soliton case above applies and the plane 
wave mode picks up the phase shift 8k associated by Levinson's theorem with the formation 
of the soliton translation modes. 
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IV. STATISTICAL PROPERTIES - CORRELATIONS - SCALING 



In this section we discuss the statistical properties of a growing interface on the basis of 
the canonical phase space formulation. Generally, we can express the noisy Burgers equation 
in Eq. (jl.ljl in the form 

du 5F 

^- = z/V— + XuVu + Vrj, (4.1) 
at ou 

where the free energy F driving the diffusive term is given by 

F = ^J dxu 2 {x). (4.2) 

For A = we have the linear Edwards- Wilkinson equation describing the temporal fluctua- 
tions in a thermodynamic equilibrium state with temperature T = A/2z//cb with stationary 
distribution given by Eq. (jl.21|) . i.e., P st oc exp[— F/ksT]. In the presence of the nonlinear 
mode coupling term XuVu Eq. (|4.1|) does not describe a thermodynamic equilibrium state 
but a stationary nonequilibrium state or kinetic growth state. It is, however, a particular 
feature of the one dimensional case that the stationary distribution is known and given 
by P st oc exp[— F /ku T], independent of A. This property also follows from the potential 



condition m 



MM 



dx[\Vu - ^^V(m 2 )] = 0, (4.3) 

which is readily satisfied since the integrand becomes a total derivative. Another way of 
noting that Eq. (|4.1|) does not describe an equilibrium state is to express the equation in 
the form du/dt = u'V5F' /Su + Vrj with effective free energy F' = (1/2) / dx(u 2 + (A/3^)w 3 ). 
Apart from the fact that F' includes an odd power in u and thus, since u = V/z, violates 
parity invariance it is also unbounded from below for u —>■ — oo and thus cannot describe a 
stable thermodynamic state. 

The stationary distribution P s t(u) is obtained in the limit t — > oo from the transition prob- 
ability P(u l — > u, t) for a pathway from the initial configuration u l to the final configuration 
u. In this limit only the linear diffusive modes for A = persist characterized by P st . This 
is consistent with the fact that the soliton contribution yields P(u l — > u,t) oc exp[— S(t)/A], 
where the action S(t) associated with the solitonic growth modes, e.g., the two-soliton con- 
figuration, typically grows linearly with t implying that the contribution to P vanishes. 
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A. Correlations in the Edwards- Wilkinson case 



In the linear case the correlation function is easily evaluated on the Langevin level from 
the stochastic Edwards- Wilkinson equation (|2.2|) and follows directly from Eq. (|2.1j) when 
averaging over the noise, see also paper II. In wavenumber-frequency space we obtain the 
Lorentzian diffusive form 

Ak 2 

(uu)(k,u)= (4.4) 
u z + [vk A y 

with diffusive poles at u — ±iuk 2 , strength A/{vk) 2 and linewidth vk 2 . We note that both 
growing, u oc exp(uk 2 t), and decaying terms, u oc exp(— uk 2 t), contribute to the stationary 
correlations; this is in accordance with the decomposition (|2.14|) . In wavenumber-time space 
we have correspondingly 

(uu)(k,t) = A e — '*"!*!, (4.5) 

and the diffusive correlations decay on a time scale determined by 1/uk 2 . For the equal-time 
correlations we obtain in particular (uu)(k,0) = A/2u, showing the spatially short ranged 
correlations in accordance with the stationary distribution (J1.21|) . For later purposes we 
also need the spectral form, see paper II, 

(uu)(x,t) = J ±±. e **e->>*n (4.6) 

In order to illustrate the method to be used later in the soliton case we evaluate here (uu) 
in the linear case on the basis of the path integral formulation Eq. (jl,14j) . Since that the 
distributions (jl.21j) and ()2.13j) factorize in wavenumber space we have in a little detail for a 
system of size L 

(uu)(k,t) oc J ndu'pdulululYlexp ^ Z^pL^I ' ) IT ex P (~XZ' M ™| 2 ) ■ ( 4 ' 7 ) 

Changing variables in P(w^,u^,t), u { n — u\ exp(— uj n t) u n , it is an easy task to carry out 
the Gaussian integrals over u n and u l n and retrieve (uu)(k,t) in Eq. (|4.5J) ; the evaluation 
of (uu) in the corresponding harmonic oscillator quantum field calculation was performed 
in paper II. Finally, evaluating Eq. ()4.6|) we infer the scaling form (jl.20|) with roughness 
exponent £ = 1/2, dynamic exponent z = 2, and scaling function, see also paper II, 

F(w) = ^[^u]- l/2 w- l/2 e- 1,iuw ,w = t/x z , (4.8) 
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defining the Edwards- Wilkinson universality class. 

Summarizing, the Edwards- Wilkinson equation describes a thermodynamic equilibrium 
state. The dynamical equilibrium fluctuations are characterized by the gapless dispersion 
law uj = vk 2 . The modes are extensive and diffusive and controlled by the characteristic 
decay time 1/uk 2 , depending on the wavenumber k. 

B. Switching and pathways in the Burgers-KPZ case 

Before we turn to the correlations in the nonlinear Burgers-KPZ case it is instructive to 
extract a couple of simple qualitative consequences of the dynamical approach. As discussed 
in Sec.|n]the propagation of a two-soliton configuration constitutes a simple growth situation 
where at each passage of the soliton pair the interface grows by a layer. Considering a pair 
configuration of amplitude u and size I the propagation velocity is v = —Am/2 and the 
associated action Si = (l/Q)v\\u\ 3 T . Across a system of size L we have \v\ = L/T, where 
T is the switching time, i.e., \u\ = 2L/XT. For the action associated with the transition 
pathway of adding a layer of thickness h = J x u(x' ,t)dt = \u\£ = 2LE/XT we then have 



We note that the thickness h does also depend on the pair-soliton size £ which does not 
enter in the action. However, the multiplicity or density of soliton pair which enters in the 
prefactor of the transition probability must depend on l/£ and we obtain qualitatively 



In the thermodynamic limit L — > oo, P — > and the switching via a two-soliton pathway is 
suppressed. At long times the action falls off as 1/T 2 . 

In the case of a switching pathway by means of two equal amplitude pair-solitons we 
obtain, correspondingly, noting that the pairs propagate with half the velocity, the action 
S2 = (1/4)51. Introducing heuristically a constant nucleation action S n uc\ associated with 
the noise-induced formation of a pair, i.e., the appropriate assignment of the noise field p, 
we have the general expression for the action associated with n pairs 



Si(T) 



Av L 3 



(4.9) 




(4.10) 



S n (T) 



1 Av L 3 



(4.11) 



2 3A 2 T 2' 



n 
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In Fig. ^1 we have plotted S versus T for n = 1 — 5 soliton pairs. Since the curves intersect 
we infer that the switching at long times takes place via a single soliton pair. At shorter 
times a switching takes place and the transition pathway proceeds by the excitation of 
multi-pair-solitons. This is clearly a finite size effect. 

A similar analysis of the soliton switching pathways in the case of the noise driven 
Ginzburg-Landau equation has recently been carried out. Here the analysis, corroborat- 
ing recent numerical optimization studies, is simpler because the soliton excitations are 
topological and have a fixed amplitude 



C. Anomalous diffusion of growth modes in the Burgers-KPZ case 

On the Langevin level the growth of the interface is a stochastic phenomena driven by 
noise. Parameterizing the growth in terms of growth modes corresponding to the propagation 
of a gas of independent pair-solitons in the slope field the dynamical approach allows a 
simple interpretation. Noting that the action associated with the pair mode is given by 
S = (i/6)z/A|w| 3 i and denoting the center of mass of the pair mode by x = (x\ + x 2 )/2 we 
have u = 2v/X = 2x/t\ and we obtain the transition probability 

P(»,«)««p(-|^), (4.12) 

for the random walk of independent pair-solitons or steps in the height profile. Comparing 
(I4.12J) with the distribution for ordinary random walk originating from the Langevin equation 
dx/dt = rj, (j]T]){f) = A5(t), P{x,t) oc exp(— x 2 /2At), we conclude that the growth mode 
performs anomalous diffusion. Assuming pairs of the same average size, the distribution 
(I4.12JI also implies the soliton mean square displacement, 

<* 2 >(t)oc(^j t 2 '\ (4.13) 

with dynamic exponent z = 3/2, identical to the dynamic exponent defining the KPZ 
universality class. This result should be contrasted with the mean square displacement 
(x 2 ) oc At 2 / Z , z = 2, for ordinary random walk. The growth modes thus perform superdif- 
fusion. This result is also obtained using the mapping of the KPZ equation to directed 
polymers in a random medium (22] . 
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The diffusion of solitons or growth modes is another signature of the stationary nonequi- 
librium state. Whereas the extended diffusive equilibrium modes for a particular wavenum- 
ber k are characterized by the stationary distribution P st oc exp(— (u/ AL)\uk\ 2 ), the ran- 
dom walk distribution of the growth modes P(x,t) oc i~ 2//3 exp( const x x 3 /t 2 ) vanishes for 
t — > oo. The growth modes or solitons disperse diffusively over the system and generate the 
stationary growth. 



D. Correlations in the Burgers-KPZ case - general 

As regard the correlations in the nonlinear Burgers-KPZ case the situation is more com- 
plex. The noisy Burgers equation (jl.lj) is not easily amenable to a direct analysis of the 
noise averaged correlations and we limit ourselves to a discussion of (uu)(x,t) within the 
canonical phase space approach. In order to evaluate the slope correlations (uu)(x,t) by 
means of Eq. (jl.l4j) . i.e., 

(uu)(x,t) = JuduW u f (x)u i (0)P(u f ,u i ,t)P st (u i ), (4.14) 

we note that the basic ingredient is the transition probability P(w f , u\t) from an initial 
configuration u 1 at time t = to a final configuration u l at time t. 



1. Sum rule 



Before continuing we observe that in the short time limit t — > it follows from the 
definition that P(u { , u 1 , t) — > 5(u { — u 1 ). The equal time correlations are thus determined by 
the stationary distribution P s t(^ 1 ) given by Eq. (jl.21j) . 



P st (w 1 ) oc exp 



— / <7.r //' (./•''- 



(4.15) 



and we have in wavenumber space (uu)(k,0) = A/2u. In wavenumber-frequency space we 
thus infer the general sum rule 



2^ 



(uu)(k, ui) 



A 
2u' 



(4.16) 



The sum rule is independent of the presence of the nonlinear growth term XuVu and thus 
is another consequence of the static fluctuation dissipation theorem which holds for the 
Burgers-KPZ equations 0, 36^ . 
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2. The transition probability 

As discussed in Sec. |H]the working hypothesis is that a growing interface at a particular 
time instant can be represented by a dilute gas of matched localized soliton excitations or 
growth modes with superimposed linear extended diffusive modes. From the analysis in 
Sec. II I II we thus have 

u(x, t) — u s (x, t) + Su(x, t), (4-17) 
p(x,t) =p.(x,t) + 6p(x,t), (4.18) 

where u s and p s (or <p s ) are given by the multi-soliton representations in Eqs. (|2.21j) and 
(l2~2~2l (or for <p B in Eq. (HOI). In the flat regions for constant slope 5u = (1/2)[6X + +SX-) 
and 5p = u5X~ are given by Eqs. (|H.5Hjl (across the soliton regions 5u and Sp vary in a 
more complicated manner as discussed in Sec. 1111(1 . 

Inserting Eqs. (|4.17(1 and ((4.18(1 in Eq. (Jl. 10(1 and using the equation of motion ()1. 6(1 
the action S decomposes in a soliton contribution S so \ and a linear contribution Sn n , S = 
Ssoi + Sim, where 

Ssoi = \ J dxdt{Vp s )\ (4.19) 
S'lin = - J dxdt{V8pf. (4.20) 

This decomposition implies that the transition probability P(u f , u 1 , t) accordingly factorizes 
like 

P{u\u\t) = P^{u\M.it)PUfo t M,t), (4-21) 

where P so \ oc exp(— S so \/ A) and P lin oc exp(— 5u n /A). Disregarding phase shift effects and 
amplitude modulations due to the dilute soliton gas, P\ in can be worked out as in the 
Edwards- Wilkinson case in Sec. HTJ yielding the expression (in wavenumber space) 

P(5u\ 5u { , t) oc exp 

with limiting distribution P st (^M f ) oc exp[— (z//A) J(dk/27i)\5u i k \ 2 } for t — > oo. 
For the soliton part we obtain inserting Eq. ((2.26(1 

(4.23) 
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v r dk \8u k — 5w 1 A .exp(— u>kt)\' 
AJ 2tt 1 - exp(-2u k t) 



(4.22) 



P(u\ w f , t) oc exp 



6A 



E 



u 



p+1 



u. 



l 9{u p 



u 



in terms of the soliton boundary values u p as depicted in Fig.HJ Note also that only the noise 
induced left hand solitons contribute to the action. We stress that the expression (|4.23jl by 
construction only holds in-between soliton collisions. In fact, at long times the appropriate 
expression for P(V, u f , t) must approach P st oc exp[— (z//A) / dxu f s (x) 2 ] in accordance with 
Eq. (|4.21j) . Likewise, the correct expression for the multi-soliton energy must vanish in the 
long time limit corresponding to the migration of the phase space orbit to the transient and 
stationary zero energy submanifolds p = and p = 2i/u, as discussed in paper III. 

3. Multi-soliton correlations - Scaling properties 

Inserting Eq. (J4.17|) in Eq. (J4.14j) the slope correlations separate in a soliton part and a 
linear (diffusive) part, 



(uu)(x,t) = (u s u s )(x,t) + (8u8u)(x,t). 



(4.24) 



Apart from phase shift and amplitude modulation effetcs due to the dilute soliton gas, the 
linear or diffusive correlations (SuSu) basically have the form given by Eq. ()4.6|) . For the 
soliton contribution (u s u s ) we obtain, inserting u s from Eq. (j2.21j) . P(u { s ,u\t) from Eq. 
(|4.23j) . and for the stationary distribution 



P st oc exp 



— ST 2 I _ I 

£S P=1 



exp 



\^\2 ^ ^ kpkp>\Xp Xpi | 
p^p' 



(4.25) 



and moreover introducing the soliton amplitude k p = (A/4z/)(w p+1 — u p ) from Eq. ()2.19|) . 



«.«.)(*, *)«(«) = E / n dkdv^M, tanh \k p \( Xp + v p t - x) tanh \k,\x, x 

P,q J I 



J|exp 



32z/ 4 
3AX' 



t\k n \ @( kn) 



n ex p 



AX 



;k n k n /(x n X n i 



(4.26) 



This formula expresses the contribution to the slope correlations from a multi-soliton con- 
figuration. It follows from the derivation that the expression only holds for times short 
compared to the soliton collision time. The initial configuration u\ at time t = propa- 
gates during time t to the the final configuration u s . The associated transition probability 
is given by Eq. (J4.23)) and the stationary distribution by Eq. (J4.25j) . The integration over 
initial and final configurations is effectuated by integrating over the amplitudes k p , the ve- 
locities v p and the soliton positions x p over a system of size L. Note that k p together with 
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-(\/2)(up+x + u p ) determine the slope u p . Likewise the dynamic partition function 



Q(t) is given by 



n(t) = J JJ dkidvidxi JJ exp 



32z/ 4 
3AX' 



't\k n \ 9( k Ti 



n ex p 



AX' 



"k n k n i(^x n x r , 



. (4.27) 



The complex form of Eqs. ()4.26|) and ()4.27|) have so far precluded a more detailed analysis. 
We can, however, in the limit of small damping extract the general scaling properties. For 
v — > we have k p — > oo and the soliton profile given by tanh |fe p |(x — v p t — x p ) converges 
to the sign function sgn(x — v p t — x p ), corresponding to a sharp shock wave. By inspection 
of Eq. ()4.26|) we note that a change of the length scale by a factor and 
be absorbed by a change of the integration variable k p , k p — ► fi 1 ^ 2 k p . In 
the action term this change of k p is finally absorbed by the scale transformation t /j, 3 / 2 . 
Consequently, for v — > we have (u s u s )(x,t) = F(t/x 3 ^ 2 ) in accordance with the general 
scaling form in Eq. (jl.2()j) . 



E. Correlations in the Burgers-KPZ case - the two-soliton sector 



In the weak noise limit A — > the action in Eq. (jl.llj) provides a general selection crite- 
rion determining the dominant dynamical configuration contributing to the distribution P. 
In the present section we propose that part of the leading growth morphology is constituted 
by a gas of two-soliton or pair configurations already analyzed in Sec. |H] In our numerical 
studies we have demonstrated that in the limit v — > the pair configuration does constitute 
a long lived quasi-particle j^. 

The evaluation of the contribution to the slope correlations from the two-soliton sector 
is straightforward and will permit a more detailed scaling analysis. Specializing the general 
expression in Eqs. ()4.26|) and ()4.27|) to the case of two solitons, i.e., a pair-soliton excitation, 
noting that k\ = —k 2 = (A/4z/)m, u 2 = u {u\ = u 3 = 0), and v\ — v 2 — v — —(\/2)u, and 
moreover considering the limit of small u, or, alternatively, using the expressions pertaining 
to the two-soliton case discussed in Sec. HH we have 

2 



(uu)(x,t)n 2 (t) = 



dudxidx 2 u 2 [sign(xi) — sign^)] x 



[signal — x 
uX 



vt) — sign(x 2 — x — vt)} x 



exp 



6A 



t\u\ 



exp 



A 



U \X 2 — X\\ 



(4.28) 
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with dynamic partition function 



n 2 (t) 



dudx\dx2 exp 



—t\u\ 

6A 1 1 



exp 



Xi\ 



(4.29) 



We note that the final configuration u is simply the initial two-soliton configuration u 1 
displaced vt along the axis without change of shape. This dynamical evolution is depicted 
in Fig. The integration over initial and final configurations is carried out by integrating 
over the soliton amplitude u, — oo < u < oo and the soliton positions x\ and x^ over a 
system of size L. The integration over the amplitude only contributes to integral when the 
pair-solitons overlap, as indicated in Fig. and we obtain by inspection of the overlap 
contribution, setting z = x — vt and £ = \x2 — x±\, the conditions x\ < z, x\ > z — £, xi < 0, 
and X\ > —£. For z > 0, i.e., x — vt > we have the overlap conditions < z < I and 
z — i < x\ < 0; for z < we obtain — i < z < and — I < x\ < z. Finally, integrating over 
the soliton position x± and the soliton pair size £ we arrive at the expression 



(uu)(x, t) 



1 / du u 2 exp 


6A 


U 


H 


exp 


V 

A 


x — vt\u 2 


£#>(«) 


L J duexp 


uX 
6A 


u\H 





(4.30) 



where the cut-off functions = / L d££exp(-(u/A)u 2 e) and cf = J L d£ exp(-(u/ A)u 2 £) 
follow from the overlap; explicitly they are given by 



C 



c 



(2), 



u 



u 



A\ 2 1_ 

J a 4 

V J u 



1-1 



1 — exp 



v 

A' 

v 

~A 



u L exp 



u 2 L 



A 



u 2 L 



(4.31) 
(4.32) 



The overall factor 1/L reflects the weight of a single pair-soliton contribution to the cor- 
relations function. In the thermodynamic limit L — > oo this contribution vanishes. For a 
dilute gas of pair-solitons of density n we expect 1/L to be replaced by n. On the other 
hand, the further L dependence of the cut-off functions, is a feature of the extended nature 
of the pair-soliton already discussed in Sec. [HJ Both and vanish as a function of 



u over a scale set by yA/uL. For u — > oo C^ 1 -* ~ 1/w 4 and ~ 1/m 2 ; for u = we have 
C« = L 2 /2 and = L. 



F. General scaling properties 



The last issue we deal with is the scaling properties of a growing interface. The dynamical 



scaling hypothesis 



17 



and general arguments based on the renormalization group fixed- 
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point structure 2J, |25[ imply the following long time - large distance form of the slope 



correlations in the stationary state: 

(uu)(x,t) = (A/2v)x 2C - 1 F(x/£(t)). (4.33) 

Here F is the scaling function and the roughness exponent ( follows from the explicitly 
known stationary distribution in Eq. (|1.21|) . the fluctuation dissipation theorem. Within 
the canonical phase space approach the stationary distribution follows from the structure 
of the zero-energy manifolds which attract the phase space orbits in the long time limit 
t — > oo, see paper III. The dynamic exponent z = 3/2 is inferred from the gapless soliton 
dispersion law in Eq. ()2.36|) . see paper II. Since the formulation is entirely Galilean invariant 
the exponent z also follows from the scaling law £ + z = 2 in Eq. (|1.3|) . 

The lateral growth of fluctuations along the interface is conveniently characterized by the 
time dependent correlation length Note that for a finite system of size L the correlation 
length saturates at the crossover or saturation time t co determined by £(£ C o) — L. In the 
linear Edwards- Wilkinson case £(t) characterizes the growth of diffusive modes and has the 
form = (z/t) 1//2 , consistent with the spectral form in Eq. (|4.6|) . In the Burgers-KPZ case 
£(£) describes the propagation of soliton modes and is given by £(£) = (A/z/) 1 / 3 (At) 2 / 3 . The 
limiting form of the scaling function lim^^oo F(w) = 1 for a; > ((t). In the dynamical regime 
for £(£) ^> x the correlation decay, i.e., (uu)(x,t) — > (u)(u) = 0, and the scaling function 
vanishes like F(w) oc w 2 ^^ for w — > 0. In Fig. ^] we have depicted the correlation length 
£(t) as a function of time for a system of size L, indicating the crossover behavior in the 
Edwards- Wilkinson and Burgers-KPZ cases. In Fig. [HI] we have plotted the time scale T as 
a function of system size indicating the various dynamic regimes. 



G. Scaling properties in the two-soliton sector 

In discussing the scaling properties associated with the two-soliton sector it is convenient 
to introduce the model parameters £q, setting the microscopic length scale, to, setting the 
microscopic time scale, t co defining the crossover or saturation time for a system of size L, 
and the correlation length £(£),note that A = £o/to, 

£ = ~, (4.34) 
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to = 4' ( 4 - 35 ) 

t co = t Q {L/i Q f'\ (4.36) 
£(t)=£o(tAo) 2/3 . (4.37) 

Rescaling the amplitude variable u we can then express the pair correlations in the form 

£ n [ duexp -fid 3 /- exp -4w 2 |f + u-^-\ F 1 (u) 
WW = i- ^— r z ,„ ; i - > ( 4 - 38 ) 



£ /ciMexp 


4 1 13^" 

Q Uf , 

O L C _ 


exp 


-4m 2 




Fi(«) 


L 


J duexp 


--Id 3 — " 

3 tea . 


F 2 (u) 



where the cut-off functions originating from the overlap are given by 

F 1 H = ^-(l + i ^)exp(-4 M 2 ), (4.39) 

F2(M) = 4^( 1_eXp(_4M2) )- (440) 

The expression (14. 38(1 holds for t > and is even in x (seen by changing u to —u). It samples 
the soliton pair propagating with velocity Am/2 and is in general agreement with spectral 
form discussed in the quantum treatment in paper II. 

The weight of a single soliton pair is of order 1/L and the correlation function (uu) thus 
vanishes in the thermodynamic limit L — ► oo. For a finite system L enters setting a length 
scale together with the saturation time t co defining a time scale, and (uu) is a function 
of x/L and t/t co as is the case for the two-soliton expression ([4.38)1 . It is instructive to 
compare this dependence with the wavenumber decomposition of (uu) in the linear diffusive 
case for A = 0. Here (uu)(xt) oc (1/L) Dn^o ex P(~ (2nn) 2 t/L 2 ) exp(innx / L) , depending on 
x/L and t/L 2 , corresponding to the saturation time t co cx L 2 , z = 2. Keeping only one 
mode for n = 1 the correlations (uu) has the same structure as in the soliton case. In the 
linear case we can, of course, sum over the totality of modes and in the thermodynamic 
limit L — » oo replace (l/L)J2 n by / dk/2n obtaining the intensive correlations (uu)(xt) = 
(A/2v)(4:iivt)~ 1 / 2 exp(— x 2 /2vt). Similarly, we expect the inclusion of multi-soliton modes 
to allow the thermodynamic limit to be carried out yielding an intensive correlation function 
in the Burgers case. 

For a finite system we have in general [(3(3] (uu)(xt) = (l/L)Gi(x/L,t/L 3 / 2 ) with 
scaling limits: Gl(x/L,0) oc const, for x ~ L, Gl(x/L,0) oc L/x lor x L and 
G L (0,t/L 3 / 2 ) oc const, for t > L 3 / 2 , G L (0,t/L 3 / 2 ) oc L/t 2 / 3 for t < L 3 ' 2 . For L -> oo 
we obtain G L (x/ 'L,t/ 'L 3 ' 2 ) -> (L/x)G(x/t 2 ' 3 ) in conformity with 
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It is an important feature of the two-soliton expression ()4.38|) that the dynamical soliton 
interpretation directly implies the correct dependence on the scaling variables x/L and 
t/t co oc t/L 3 ^ 2 , independent of a renormalization group argument. However, the scaling 
limits are at variance with Gl- Setting, according to ()4.38|) (uu)(xt) = (£ /L)F(x/L,t/t co ), 
F(x/L, 0) assumes the value .47 for x -C L and decreases monotonically to the value ~ .08 
for x ~ L, whereas Gl diverges as L/x for x<L. Likewise, F(0,t/t co ) decays from .47 for 
t <C t co oc L 3 / 2 to for t ^> t co ; for t ~ t co we have F 2 ~ .15, whereas Gl diverges as L/i 2 / 3 
for t <^:t co . 

This discrepancy from the scaling limits is a feature of the two-soliton contribution which 
only samples the correlation from a single soliton pair. Moreover, at long times the soliton 
contribution vanishes and the scaling function is determined by the diffusive mode con- 
tribution in accordance with the convergence of the phase space orbits to the stationary 
zero-energy manifold. We note, however, the general trend towards a divergence for small 
values of x and t is a feature of F. 

Introducing the scaling variables w = x/£ oc xjt 2 !' 3 and r = t/t s oc t/L 3 ^ 2 we can also 
express (14.38)1 in the form 



(uu)(xt) = (£ /L)F(w,r) , (4.41) 
where the scaling function F is now given by 



u 



F(w, t) = J - j— ^ , (4.42) 

fdue-iM* T F 2 (u) 

and summarize our findings in Fig. where we have depicted F(w,r) for a range of r 
values. For fixed small w = x/£ oc xjt 2 !' 3 we have F — > .47 for r = t/t s oc t/L 3 ^ 2 — > 0; 
for large r we obtain F — > 0. The weak maximum moving towards smaller values of u> for 
decreasing r is a feature of the functional form of F in ()4.42j) and thus due to the soliton 
approximation. The true scaling function is not expected to have any particular distinct 
"3, 27, 28, 3, 



features 



22 



24, 



30, 31 
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V. SUMMARY AND CONCLUSION 

In the present paper we have continued our analysis of the noisy Burgers equation in 
one spatial dimension within the weak noise canonical phase space approach developed in 
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previous papers. We believe that the noisy Burgers equation or the equivalent KPZ equation, 
which have been studied intensively, is of fundamental and paradigmatic significance in the 
context of a continuum field theoretical description of nonlinear non equilibrium phenomena. 
The advantage of the canonical phase space method which is an elaboration and a dynamical 
system theory interpretation of the saddle point equations originating from the Martin- 
Siggia-Rose functional formulation or, equivalently, a phase space formulation of the Freilin- 
Wentzel variational approach to the Fokker-Planck equation, actually dating back to work 



by Machlup and Onsager |61|,|62|, is that it replaces the stochastic Langevin equation with 
coupled deterministic field equations, yielding on the one hand an interpretation of the 
growth morphology and pattern formation and on the other hand a practical scheme for the 
evaluation of the statistical properties and correlations in the weak noise limit. 

Here we have in some detail discussed i) the growth morphology engendered by the 
propagation of domain walls or solitons, the growth modes, ii) the superimposed linear 
modes and their transmutation to propagating modes in the presence of the growth modes, 
and, finally, iii) the statical and scaling properties, particularly, in the two-soliton sector. 
The weak noise theory of the one dimensional Burgers or KPZ equation is, however, far from 
complete and many open questions remain. We mention below a series of topics which it 
would be of considerable interest to investigate: i) the interpretation of the solitoninc growth 
picture in the context of the mapping of the KPZ equation to the model of directed polymers 
in a random medium, ii) a more complete analysis of the multi-soliton correlations in the 



thermodynamic limit with the purpose of making contact with ot 
universality class, e.g. the polynuclear growth model 



ier models in the KPZ 
iii) elaboration of the 



anomalous diffusion of growth modes, iv) contact with other models for many-body systems 
far from equilibrium, e.g. driven lattice gas models 5^|, and, finally, iv) the extension of the 
weak noise approach to higher dimensions. 
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FIG. 1: We depict the noisy behavior of a wavenumber component of the slope field, Uk(t), in 
the Edwards- Wilkinson case for A = 0. After a transient period given by 1/cjfc the noise on the 
same time scale gradually picks up the motion and drives Uk{t) into a stationary noisy state. The 
transient regime is denoted by I; the stationary regime by II. 



FIG. 2: Canonical phase space plot in the linear case for A = of a single wavenumber component. 
The solid lines indicate the transient submanifold p^ = (I) and the stationary submanifold 
Pk = 2uuk (II). The stationary saddle point (SP) is at the origin. For t — ► oo the orbit from u\ to 
u\. migrates to the zero-energy manifold. The infinite waiting time at the saddle point corresponds 
to ergodic behavior. 
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FIG. 3: We depict the right hand and left hand moving solitons forming the "quarks" in the 
description of a growing interface. We have, moreover, shown the associated height profiles. 
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FIG. 4: We depict an n-soliton slope configuration of a growing interface. The p-th soliton moves 
with velocity v p = — (A/2)(u p +i +u p ), has boundary values u + and u_, and is centered at x p . The 
arrows on the horizontal inter-soliton segments indicate the propagation of linear modes. 
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FIG. 5: We depict the 4-soliton representation of the slope field u, the noise field p, and the height 
field h. The shaded area in u, i.e., the integration of u up to the point x equals the height h at x. 
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FIG. 6: We depict the slope field 112, the associated noise field P2, and the resulting height profile 
I12 at time t = for a two-soliton configuration. The slope configuration has amplitude u, size 
I = \x2 — xi\, and propagates with velocity v = —An/2. The arrow indicates the propagation of 
the superimposed linear mode with phase velocity 2v (discussed in Sec. IHIj) . 
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FIG. 7: The complex eigenvalue spectrum for the fluctuating linear modes, ReEk, K = u(k 2 + kg — K 2 ) 
and lmEk ;K = 2vkn. The bounding parabola for k = k s and k = —k s corresponds to the left and 
right extended modes propagating towards the soliton center; they are denoted I and I', respectively. 
The shaded area bounded by the parabola corresponds to localized propagating modes for 
For k = the spectrum is real corresponding to a localized non-propagating symmetric mode. The 
point TM corresponds to the translation mode. 



55 





8u 






FIG. 8: The fluctuation patterns of the pinned dynamical modes corresponding to the sectors of 
the eigenvalue spectrum in Fig. [7| The arrows indicate the propagation directions. 
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FIG. 9: We depict the extended mode propagating from right and left towards the soliton center 
which acts like a sink. The center point at x = oscillates with frequency uj = kv. The arrows 
indicate the propagation direction. 
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FIG. 10: The action given by Eq. (|4.11|) is plotted as a function of T for transition pathways 
involving up to n = 5 soliton pairs. The lowest action and thus the most probable transition 
is associated with an increasing number of soliton pairs at shorter times, indicated by the heavy 
limiting curve. The curves are plotted in arbitrary units. 
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FIG. 11: We depict the two-soliton configuration in the limit v — > contributing to the slope 
correlations (uu). The initial pair u 1 propagates to the final configuration in time t with velocity 
v = -An/2. 
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FIG. 12: The two-soliton configuration of size I = \xi — X2\ and amplitude u. The shaded overlap 
area of size U — x yields a contribution to the slope correlation function. 
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FIG. 13: The correlation lengths £ew(*) = {ut) 1 / 2 and f B (t) = (A//^) 1/3 (At) 2 / 3 as functions of t. 
For a finite system of size L the correlation lengths define the crossover times i^ff oc L 2 /u and 
t^o 1 fx A _1 (^/A) 1 / 2 L 3 / 2 , determining the transition from transient to stationary growth. 
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FIG. 14: In the early time regime for T <C t s c ^ the distribution is dominated by solitons. In the 
intermediate time regime for t^ ff » T 2> i™ 1 the solitons become suppressed and are replaced by 
the diffusive modes. Finally, for T » t^ ff the diffusive modes also die out and we approach the 
stationary distribution. 
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FIG. 15: Plot of the scaling function F(w, r) as a function of the scaling variable w = x/£ oc x 

/t 2/3 

for a range of values of r = t/t co oc i/L 3 / 2 . 



59 



